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We analyze the cosmological implications of F-term hybrid inflation with a subdominant Fayet- 
Iliopoulos D-term whose presence explicitly breaks a D-parity in the inflaton-waterfall sector. 



Q^' This scenario of inflation, which is called -Fo-term hybrid model for brevity, can naturally predict 
^ I lepton number violation at the electroweak scale, by tying the /i-parameter of the MSSM to an 
S0(3)-symmetric Majorana mass m^v, via the vacuum expectation value of the inflaton fleld. We 
^ i show how a negative Hubble-induced mass term in a next-to-minimal extension of supergravity 
helps to accommodate the present CMB data and considerably weaken the strict constraints 
on the theoretical parameters, resulting from cosmic string effects on the power spectrum P-ji. 
The usual gravitino overabundance constraint may be signiflcantly relaxed in this model, once 
the enormous entropy release from the late decays of the ultraheavy waterfall gauge particles 
is properly considered. As the Universe enters a second thermalization phase involving a very 
low reheat temperature, which might be as low as about 0.3 TeV, thermal electroweak-scale 
resonant leptogenesis provides a viable mechanism for successful baryogenesis, while thermal right- 
handed sneutrinos emerge as new possible candidates for solving the cold dark matter problem. 
In addition, we discuss grand unifled theory realizations of i</)-term hybrid inflation devoid of 
cosmic strings and monopoles, based on the complete breaking of an SU(2)x subgroup. The 
i*£)-term hybrid model offers rich particle-physics phenomenology, which could be probed at high- 
energy colliders, as well as in low-energy experiments of lepton flavour or number violation. 
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1 Introduction 



Standard big-bang cosmology faces severe difficulties in accounting for the observed flatness 
and enormity of the causal horizon of today's Universe. It also leaves unexplained the 
origin of the nearly scale-invariant cosmic microwave background (CMB), as was found by 
a number of observations over the last decade [1-5]. All these pressing problems can be 
successfully addressed within the field-theoretic framework of infiation [6] . As a source of 
infiation, it is usually considered to be a scalar field, the infiaton, which is displaced from 
its minimum and whose slow-roll dynamics leads to an accelerated expansion of the early 
Universe. In this phase of accelerated expansion or inflation, the quantum fluctuations of 
the inflaton fleld are stretched on large scales and eventually get frozen when they become 
much bigger than the Hubble radius. These quantum fluctuations get imprinted in the 
form of density perturbations, when the former are crossing back inside the Hubble radius 
long after inflation has ended. In this way, inflation provides a causal mechanism to explain 
the observed nearly-scale invariant CMB spectrum. 

A complete description of the CMB spectrum involves about a dozen of cosmological 
parameters, such as the power spectrum P-ji of curvature perturbations, the spectral index 
^s, the running spectral index dus/dlnk, the ratio r of tensor-to-scalar perturbations, the 
baryon-to- photon ratio of number densities tje, the fractions of relic abundance Qdm and 
dark energy Q\ and a few others. Recent WMAP data [2,4], along with other astronomical 
observations [3] , have improved upon the precision of almost all of the above cosmological 
observables. In particular, the precise values of these cosmological observables set stringent 
constraints on the model-building of successful models of inflation. To ensure the slow-roll 
dynamics of the inflaton, for example, one would need a scalar potential, which is almost 
flat. Moreover, one has to assure that the flatness of the inflaton potential does not get 
spoiled by large quantum corrections that depend quadratically on the cut-off of the theory. 
In this context, supersymmetry (SUSY), softly broken at the TcV scale, emerges almost as 
a compelling ingredient not only in the model-building of inflationary scenarios, but also 
for addressing technically the so-called gauge-hierarchy problem. 

One of the most predictive and potentially testable scenarios of inflation is the model 
of hybrid inflation [7] . An advantageous feature of this model is that the inflaton may start 
its slow-roll from fleld values well below the reduced Planck mass m-p\ — 2.4 x 10^® GeV. 
As a consequence, cosmological observables, such as P-ji and Ug, do not generically receive 
significant contributions from possible higher-dimensional non-renormalizable operators, 
as these are suppressed by inverse powers of 1/mpi. Thus, the hybrid model becomes 
very predictive and possibly testable, in the sense that the infiaton dynamics is mainly 
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governed by a few renormalizable operators which might have observable imphcations for 
laboratory experiments. In the hybrid model, inflation terminates through the so-called 
waterfall mechanism. This mechanism is triggered, when the infiaton field passes below 
some critical value 0c- In this case, another field X different from 0, which is called the 
waterfall field and is held fixed at origin initially, develops a tachyonic instabifity and rolls 
rapidly down to its true vacuum expectation value (VEV). 

Hybrid inflation can be realized in supersymmetric theories in two forms. In the first 
form, the hybrid potential results from the F-terms of a superpotential, where the slope 
of the potential may come either from supergravity (SUGRA) corrections [8] and/or from 
radiative effects [9]. The second supersymmetric realization [10] of hybrid inflation uses a 
dominant Fayet-Iliopoulos (FI) D-term [11], which may originate from an anomalous local 
U(1)q symmetry within the context of string theories. 

All models of inflation embedded in SUGRA have to address a serious problem. 
This is the so-called gravitino overabundance problem. If abundantly produced in the 
early Universe, gravitinos may disrupt, via their late gravitationally-mediated decays, the 
nucleosynthesis of the light elements. In order to prevent this from happening, gravitinos 
G must have a rather low abundance today, i.e. = n^/s < lO'-'^^-lO^"'^^, where 
is the number density of gravitinos and s is the entropy density. The upper bound on 
Yg depends on the properties of the gravitino and becomes tighter, if gravitinos decay 
appreciably to hadronic modes. These considerations set a strict upper bound on Universe's 
reheat temperature Treh, generically implying that Treh ~ 10^°-10^ GeV [51,56]. This upper 
limit on Treh severely restricts the size of any renormalizable superpotential coupling of 
the infiaton to particles of the Standard Model (SM). All these coupfings must be rather 
suppressed. Typically, they have to be smaller than about 10~^ [12]. 

The aforementioned gravitino constraint may be considerably relaxed, if there is a 
mechanism that could cause late entropy release in the evolution of the early Universe. Such 
a mechanism could then dilute the gravitinos to a level that would not upset the limits 
derived from Big Bang nucleosynthesis (BBN). This possibility might arise even within 
the context of F-term hybrid infiation, if a subdominant FI D-tadpole associated with the 
gauge group U(l)x of the waterfall sector were added to the model. Such a scenario was 
recently discussed in [13]. It has been observed that the presence of a FI D-term breaks 
explicitly an exact discrete symmetry acting on the gauged waterfall sector, i.e. a kind 
of D-parity, which would have remained otherwise unbroken even after the spontaneous 
symmetry breaking (SSB) of U(l)x- As a consequence, the ultraheavy U(l)x-gauge-sector 
bosons and fermions, which would have been otherwise stable, can now decay with rates 
controlled by the size of the FI D-term. Since these particles could be abundantly produced 
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during the preheating epoch, their late decays could give rise to a second reheat phase in 
the evolution of the early Universe. Depending on the actual size of the FI Z^-term, this 
second reheat temperature may be as low as 0.5-1 TeV, resulting in an enormous entropy 
release. This could be sufficient to render the gravitinos under abundant, which might be 
copiously produced during the first reheating from the perturbative infiaton decays. 

In this paper we present a detailed analysis of F-term hybrid inflation with a sub- 
dominant FI D-tadpole. As mentioned above, the presence of the FI Z^-tadpole is essential 
for explicitly breaking an exact discrete symmetry, a D-parity, which was acting on the 
gauged waterfall sector. In [13], we termed this inflationary scenario, in short, Fp-term 
hybrid inflation. As the inflaton chiral superfleld S couples to the Higgs-doublet chiral 
superfields Hu,d, through XSHuHd, the model generates an effective /^-parameter for the 
Minimal Supersymmetric Standard Model (MSSM), through the VEV {S) [14]. The same 
mechanism may also generate an effective Majorana mass matrix for the singlet neutrino 
superfields A^i,2,3 [15], through the operator ^ pij SNiNj [13,16]. Assuming that this last 
operator is S0(3)-symmetric or very close to it, i.e. pij ~ p I3, the resulting lepton-number- 
violating Majorana mass, niN — p (S), will be closely tied to the ^u-parameter of the MSSM. 
If A ~ p, the Fo-term hybrid model will then give rise to 3 nearly degenerate heavy 
Majorana neutrinos J^i,2,3ii, as well as to 3 complex right-handed sneutrinos A^i,2,3, with 
electroweak-scale masses. Such a mass spectrum opens up the possibility to explain the 
baryon asymmetry in the Universe (BAU) r)B [17,18] by thermal electroweak-scale resonant 
leptogenesis [16,19,20], almost independently of the initial baryon- number composition of 
the primordial plasma. Moreover, since the F^-term hybrid model conserves it!-parity [13], 
the hghtest supersymmetric particle (LSP) is stable. Here, we examine the possibihty that 
thermal right-handed sneutrinos are responsible for solving the cold dark matter (CDM) 
problem of the Universe. 

In this paper we also improve an earlier approach [13], concerning the production of 
the quasi-stable U(l)x gauge-sector particles during the preheating epoch. In addition, 
we present a numerical analysis that properly takes into account the combined effect on 
the reheat temperature Tj-eh from the infiaton and gauge-sector particle decays and their 
annihilations. We call this two-states' mechanism of reheating the Universe, coupled reheat- 
ing. After solving numerically a network of Boltzmann equations (EEs) that appropriately 
treat coupled reheating, we obtain estimates for the present abundance of gravitinos in the 
Universe. We show explicitly, how a small breaking of D-parity sourced by a subdominant 
FI £)-tadpole helps to relax the strict gravitino overproduction constraint. 

In addition to gravitinos, one might have to worry that topologically stable cosmic 
strings do not contribute significantly to the CMB power spectrum P-jz. Cosmic strings. 
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global or local, are topological defects and usually form after the SSB of some global or local 
U(l) symmetry [21-23]. According to recent analyses [24], cosmic strings, if any, should 
make up no more than about 10% of the power spectrum Pji. This last requirement puts 
severe hmits on the allowed parameter space of models of inflation. There have already 
been some suggestions [25] on how to get rid of cosmic strings, based on modified versions 
of hybrid infiation. Here, we follow a different approach to solving this problem. We 
consider models, for which the waterfall sector possesses an SU(2)x gauge symmetry which 
breaks completely, i.e. SU(2)x — > I, such that neither cosmic strings nor monopoles are 
produced at the end of infiation. In this case, gauge invariance forbids the existence of 
an SU(2)x D-tadpole D"^. However, Planck-mass suppressed non-renormalizable operators 
that originate from the superpotential or the Kahler potential can give rise to explicit 
breaking of £)-parity. The latter may manifest itself by the generation of effective D"^- 
tadpole terms that arise after the SSB of SU(2)x . In this way, all the SU(2)x gauge-sector 
particles can be made unstable. 

The organization of the paper is as follows: in Section 2, we describe the Fo-term 
hybrid model and calculate the 1-loop effective potential relevant to infiation. In addition, 
we discuss the possible cosmological consequences of radiative effects on the fiat directions 
in the MSSM. We conclude this section by outhning how the Fp-term hybrid model could 
generally be embedded into a grand unified theory (GUT), including possible reahzations of 
a GUT without cosmic strings and monopoles. Technical details concerning mechanisms of 
explicit D-parity breaking in SUGRA, e.g. via an effective subdominant D-tadpole or non- 
renormalizable operators in Kahler potential, are given in Appendix A. Section 3 analyzes 
the constraints on the theoretical parameters, which are mainly derived from considerations 
of the power spectrum P-jz and a strongly red-tilted spectral index Us, with Us ~ 0.95, as 
observed most recently by WMAP [4,5]. We show how a negative Hubble-induced mass 
term in a next-to-minimal extension of supergravity helps to account for the present CMB 
data, as well as to substantially weaken the strict constraints on the model parameters, 
originating from cosmic string effects on P-jz, within a U(l)x realization of the Fo-term 
hybrid model. 

In Section 4, we analyze the mass spectrum of the infiaton- waterfall sector in the post- 
infiationary era and present naive estimates of the reheat temperature T^eh as obtained from 
perturbative infiaton decays. We then make use of an improved approach to preheating 
and compute the energy density of the quasi-stable waterfall gauge particles. In Section 5, 
we solve numerically the BEs relevant to coupled reheating and present estimates for the 
gravitino abundance in the present Universe. In Section 6, we demonstrate, how thermal 
electroweak-scale resonant leptogenesis can be realized within the Fp-term hybrid model 
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and discuss the possibility of solving the CDM problem, if thermal right-handed sneutrinos 
are considered to be the LSPs in the spectrum. In Section 7, we present our conclusions, 
including a summary of possible particle-physics implications of the F^-term hybrid model 
for high-energy colliders and for low-energy experiments of lepton flavour and/or number 
violation. 



2 General Setup 

In this section we first present the general setup of the Fp-term hybrid model within 
the minimal SUGRA framework and compute the renormalized 1-loop effective potential 
relevant to inflation. We then discuss the cosmological implications of radiative effects on 
the MSSM flat directions for F^-term hybrid inflation and for SUSY inflationary models 
in general. Finally, we analyze the prospects of embedding the F^-term hybrid model into 
a GUT. 



2.1 The Model 

The renormalizable superpotential of the F^-term hybrid model is given by 

W = kS {XiX2 - M^) + XSH^Hd + ^SNiNj + h'^^L^H^Nj 

+ W^mTsS , (2.1) 

where W]^ssm denotes the MSSM superpotential without the /x-term: 

^mTs°m = hlQiHjj^ + h%HaQiD^ + hiHoLiEi . (2.2) 

The first term in (2.1) describes the infiaton-waterfaU (IW) sector. Specifically, S is the 
SM-singlet infiaton superfield, and Xi^2 is a chiral multiplet pair of the waterfall fields with 
opposite charges under the U(l)x gauge group, i.e. Q{Xi) — — (5(^2) = 1- In addition, 
the corresponding infiationary soft SUSY-breaking sector obtained from (2.1) reads: 

-£soft = MlS*S+(KA,SX^X2 + \AxSH^Ha + ^A,SN,Ni-KasM^S + (2.3) 

where M5, ^k,a,p and as are soft SUSY-breaking mass parameters of order Msusy ~ 1 TeV. 

The second term in (2.1), A SHuHd, induces an effective /^-parameter, when the scalar 
component of S, S, acquires a VEV, i.e. 

Pi ^ X{S) ^ ^\A,-as\ . (2.4) 
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In obtaining the last approximate equality in (2.4), we neglected the VEVs of Hu^d and 
considered the fact that the VEVs of the waterfall fields X1.2 after inflation are: {Xi^2) = 
M [14]. For A ~ the size of /x-parameter turns out to be of the order of the soft- 
SUSY breaking scale Msusy, as required for a successful electroweak Higgs mechanism. 
By analogy, the third term in (2.1), ^ pij S NiNj, gives rise to an effective lepton-number- 
violating Majorana mass matrix, i.e. Ms = pij vs- Assuming that pij is approximately 
S0(3) symmetric, viz. pij fsi p I3, one obtains 3 nearly degenerate right-handed neutrinos 
^1,2,3 R, with mass 

rriN = pvs ■ (2.5) 

If A and p are comparable in magnitude, then the //-parameter and the S0(3)-symmetric 
Majorana mass itln are tied together, i.e. mv ~ /W, thus leading to a scenario where the 
singlet neutrinos 1^1^2,3 r can naturally have TcV or electroweak-scale masses [13,16]. 

The renormalizable superpotential (2.1) of the model may be uniquely determined by 
imposing the continuous R symmetry: 

S ^ e'^'S, L ^ e'^^L, Q ^ e''' Q , (2.6) 

with W e^°'W, whereas all other fields remain invariant under an R transformation. 
Notice that the R symmetry (2.6) forbids the presence of higher- dimensional operators of 
the form XiX2NiNj/mpi. This fact ensures that the electroweak-scale Majorana mass mjv 
does not get destabilized by Planck-scale SUGRA effects. 

One may now observe that the superpotential (2.1) is symmetric under the permuta- 
tion of the waterfall fields, i.e. Xi <-> X2. This permutation symmetry persists, even after 
the SSB of U(l)x, since the ground state, (Xi) = {X2) — M . is invariant under the same 
symmetry as well. Hence, there is an exact discrete symmetry acting on the gauged water- 
fall sector, a kind of D-parity. As a consequence of /^-parity conservation, the ultraheavy 
particles of mass gM, which are related to the U(l)x gauge sector, are stable. Such a 
possibility is not very desirable, as these particles, if abundantly produced, may overdose 
the Universe at late times. In order to break this unwanted /^-parity, a subdominant FI 
D-term, —^grripiD, is added to the model [13], giving rise to the D-term potential ^ 

Vd = J {\X,f - |X2|2 - ml,y . (2.7) 

-'^The D-parity is an accidental discrete symmetry and it should not be confused with the U(l)j^ charge 
conjugation symmetry realized by the transformations: Xi and X2 ^ Although both discrete 

symmetries have the same effect when acting on the U(l)j^ scalar current = K-^i '9^ ^1 ~ -^2 ^2), 
i.e. ^ ~jx^ t^^y crucially differ when they are apphed on the FI D-term: — |to|j(|XiP — |X2|^). This 
term is even under charge conjugation, but odd under a £>-parity conjugation. 
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The FI D-tcrm will not affect the inffationary dynamics, as long as gmpi <^ kM. Techni- 
cally, a subdominant D-term can be generated radiatively after integrating out Planck-scale 
heavy degrees of freedom. Further discussion is given in Section 4.1 and in Appendix A, 
where we also discuss the possibility of breaking expUcitly /^-parity by non-renormalizable 
Kahler potential terms. The post- inflationary imphcations of the FI £)-term, mpi, for the 
reheat temperature Treh and the gravitino abundance Yq will be analyzed in Section 5. 

The inflationary potential V^inf may be represented by the sum 



^nf = l^i? + -^i? + I^SUGRA , 



(2.8) 



where V-^^^ and V-^l^ are the tree-level potential and the 1-loop effective potential, respec- 
tively and VsuGRA contains the SUGRA contribution. Including soft-SUSY breaking terms 
related to S, the tree-level contribution to the inflationary potential is 



ZsK^M^ + MlS*S - {nasM'^S + H.c.) , 



(2.9) 



1/2 

where Zg is the wave-function renormalization of the inflaton field which is needed to 
renormalize the 1-loop effective potential given below in the SUSY limit of the theory. The 
counter-term, 5Zs = Zs — 1, due to S wave-function renormalization may be obtained from 
the inflaton self-energy Ussip'^), through the relation 



SZc = 



dReUssip') 



(2.10) 



Calculating the UV part of SZs from this very last relation, we find 



SZc 



1 



327r2 



2J\fK^ In 



+ In 



+ In 



(2.11) 



= M is taken 



where is the renormalization scale and the infiaton field value 15*^1 
common mass renormalization point. In addition, the parameter A/" in (2.11) represents 
the dimensionality of the waterfall sector. For example, it is A/" = 1 for an U(l)x waterfall 
sector, whilst it is A/" = A^, if Xi (X2) belongs to the fundamental (anti-fundamental) 
representation of an SU(A^) theory. Observe, finally, that only the fermionic components 
of the superfields, Xi 2, Hu,d, A^i,2,3, contribute to 5Zs. 

Ignoring soft SUSY-breaking terms, the 1-loop effective potential relevant to infiation 
is calculated to be 



327r2 
+ 2A^ 



\S'' + M''\^ In 



2|2 



In 



A2(|5|2 + fM2) 



+ S' - M'\' In 



+ - f M' 



2|2 



In 



A2(|5|2-£M2; 
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\S^ M^r In — + 5*^ - ^M^r In ' ^ 



Given (2.11) and (2.12), it can be checked that the expression V^^£^ + V-^^^ is independent 
of In (5^, as it should be. 

Finally, the SUGRA contribution Vsugra to l^nf in (2.8) is highly model-dependent. 
In general, one expects an infinite series of non-renormalizable operators to occur in T^ugra, 
i.e. [8,26,27] 



^suGRA = -c'hH'\S\' + k'M'^ + 0{\S\') . (2.13) 

2 mpj 



where = /t^M^/(3mpj) is the squared Hubble rate during inflation. The flrst term 
in (2.13) represents a Hubblc-induccd mass term, which is preferably deflned to be negative 
for observational reasons to be discussed in Section 3. In a model with a minimal Kahler 
potential, the parameter ch vanishes identically.^ In fact, if \ch\ ~ 10~^, its influence on 
the CMB data [29] gets marginalized. In our analysis in Section 3, we present results for 
two representative models: (i) the scenario with a minimal Kahler potential {ch = 0); (ii) a 
next-to-minimal Kahler potential scenario with Ch ~ 0.2, where only the effect of the term 
{S'^Sy/mp^ is considered and all higher order non-renormalizable operators are ignored 
in the Kahler manifold. Moreover, we neglect possible 1-loop contributions to Vinf from 
^K,A,p-terms, which are insignificant for values M > 10^^ GeV. We only include the tadpole 
term nasM'^ S, which may become relevant for values of k < 10~^, but ignore all other soft 
SUSY-breaking terms, since they are negligible during infiation [12]. 

The stability of the inflationary trajectory in the presence of the Higgs doublets Hu^d 
and the right-handed scalar neutrinos provides further restrictions on the couplings 

A and p. In order to successfully trigger hybrid inflation, the flelds at the start of inflation 
should obey the following conditions: 

Res'- = > M, = 0, //S, = 0, Ni^,^, = 0. (2.14) 

The precise start values of the inflaton ReS*™ are determined by the number of e-folds 
Me, which is a measure of Universe's expansion during inflation (see also our discussion in 



^Strictly speaking, curvature effects related to an expanding de Sitter background will contribute to 
the potential a term given by {2Uk'^ + + 3^^) H'^\S\'^ \n{\S\'^ /Q'^), even in the minimal Kahler 

potential case [28]. Such a term, however, turns out to be negligible to affect the inflation dynamics in the 
Fo-term hybrid model. Finally, this term may be partially absorbed into the RG running of (?jj{Q'^). 
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Section 3). After inflation and the waterfall transition mechanism have been completed, it 
is important to ensure that the waterfall fields acquire a high VEV, i.e. X^''^^ = M, while 
all other fields have small electroweak-scale VEVs. This can be achieved by requiring that 
the Higgs-doublet and the sneutrino mass matrices stay positive definite throughout the 
inflationary trajectory up to a critical value \Sc\ ~ M. Instead, the corresponding mass 
matrix of Xi^2 will be the first to develop a negative eigenvalue and tachyonic instability 
close to \Sc\. As a consequence, the fields Xi^2 will be the first to start moving away from 
and set in to the 'good' vacuum Xf^^ — = M, well before the other fields, e.g. 

and 7Vi° ,3, go to a 'bad' vacuum where X^f = 0, H^f = f M and TVj" ,3 = ^ M. To 
better understand this point, let us write down the mass matrix in the weak field basis 
{Hd, H*): 

^""'^^ ~ [-kX{m'-x*x*) x'\s\' J ■ ^^-^^^ 

Then, positive definiteness of -/V^mggs implies that 

A|5|2 > k\M^ -XiX2\ . (2.16) 

Prom (2.16), it is evident that the condition X> nis sufficient for ending hybrid infiation to 
the 'good' vacuum. Finally, one obtains a condition analogous to (2.16) from the sneutrino 
mass matrix, which is equivalent to having p ^ k. The above two constraints on A and p, 
i.e. X, p > K, will be imposed in the analysis presented in Section 3. 



2.2 Radiative Lifting of MSSM Flat Directions 

Flat directions in supersymmetric theories, e.g. in the MSSM [30], play an important role 
in cosmology [31,32]. As we will demonstrate in this section, however, their infiuence on 
Fo-term hybrid infiation is minimal under rather realistic assumptions. 

One possible consequence of fiat directions could be the generation of a primordial 
baryon asymmetry 77^ through the Affleck-Dine mechanism [31]. However, if this initial 
baryon asymmetry 77^ is generated at temperatures T > mjv, it will rapidly be erased by 
the strong (B — L)-violating interactions mediated by electroweak-scale heavy Majorana 
neutrinos at T ~ niN- The BAU will then reach the present observed value by means of 
the thermal resonant leptogenesis mechanism and will only depend on the basic theoretical 
parameters of the Fo-term hybrid model [16,20]. More details are given in Section 6. 

In addition, one might argue that large VEVs associated with quasi-flat directions 
in the MSSM would make all MSSM particles so heavy after inflation, such that all per- 
turbative decays of the inflaton would be kinematically blocked and hence the Universe 
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would never thermalize [33]. The system may fall into a false vacuum with a large VEV at 
the start of inflation, which could, for example, be triggered by a negative Hubble-induced 
squared mass term of order [34] , along the flat direction. In the F^-term hybrid model, 
however, spontaneous SUSY breaking due to a non-zero {S) is communicated radiatively to 
the MSSM sector, via the renormalizable operators XSHuHd and pSNiNi. Consequently, 
their effects on the MSSM flat directions can be large and so affect the inflaton decays which 
proceed via the same renormalizable operators. In the following, we will present a careful 
treatment of this radiative hfting of MSSM flat directions, and examine the conditions, 
under which the directions would remain sufficiently flat so as to prohibit the Universe 
from thermal equilibration, shortly after inflation. 

To obtain a flat direction in supersymmetric theories, one has to impose the conditions 
of D- and F-flatness on the scalar potential V, namely the vanishing of all F- and D-terms 
for a speciflc fleld conflguration a. D-flatness is automatic for any flat direction associated 
with a gauge-invariant operator, which is absent in the MSSM, e.g. DiDjUk- Based on this 
observation, let us therefore consider here the gauge-covariant fleld conflguration 

V3 V \^Rk\ \URk\ J 

where i ^ j- It can be straightforwardly checked that the fleld conflguration o", with the 
constraint 

^ d*ni - ^ 4, - unk y^O (2.18) 

\URk\ \URk\ 

and all remaining fields being set to zero, is a flat direction, with vanishing F- and D- 
terms. It is then easy to verify that the scalar potential V{a) is truly flat, i.e. dV/da — 0. 
Although we will consider here the case of a — URk, the discussion of other squark and 
slepton flat directions is completely analogous. For notational convenience, we drop all 
generation indices from the flelds, and denote the flat direction simply by ur. 

Because of the spontaneous SUSY breaking induced by the non-zero VEV of S*, the 
flatness of the potential along the Wij-direction gets lifted, once radiative corrections are 
taken into account. The non-renormalization theorem related to theories of SUSY is still 
applicable and entails that this radiative lifting should be UV flnite and therefore calculable. 
We start our calculation by considering the pertinent mass spectrum in the background 
of a non-zero S and ur. The fermionic sector consists of 2 Dirac higgsino doublets, with 
squared masses m| — \^\S\'^ -\- h'^\uR\'^ , while the mass spectrum of the bosonic sector may 
be deduced by the mass matrix 

/ X^\S\^ -kXM^ hXSu*ji \ 

-kXM^ X^\S\'^ + h'^\uR\'^ , (2.19) 

^ hXSuR /i^lwRp / 



Ml 
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which is defined in the weak basis {Hd-, H*, Q). The couphng h in (2.19) represents a 
generic up-type quark Yukawa couphng. 

In the renormahzation scheme of dimensional reduction with minimal subtraction DR [35], 
the 1-loop efi^ective potential V^^^ related to ur is given by 

'M^\ 3" 



In 



(2.20) 



where STr denotes the usual supertrace, e.g. STrAl^ = Tr Mjj - 2m|, STrM"^ = TrMjj - 
2m4 etc. In the absence of soft SUSY-breaking terms, one finds that STrAl^ = and 
STrAl^ = 2«;^A^M^. The first condition implies the absence of quadratic UV divergences 
in SUSY theories, whereas the first together with the second one ensure the UV finiteness 
along the ur direction, namely the fact that dV^^\uR)/duR is independent. 

It would be more illuminating to compute the 1-loop effective potential in (2.20) in a 
Taylor series expansion with respect to h'^\uR\'^. To order h'^\uR\'^, the 3 mass eigenvalues 
of A1|^ are approximately given by 



Notice that in the limit ur 0, one obtains: m| = A^jS'p, M| = X\S\'^ ± kXM"^ and 
Mq = 0, as expected. Moreover, it is not difficult to check that STrTVl^ = 0{h^\uR\^) and 
STrTW^ = 2/t^A^M^ + (9(/i®|-U/j|^), in accordance with our discussion given above. 

Employing the fact that \S\'^ ^ ^ at the start of infiation, the 1-loop effective 
potential V^^\ur) may further be approximated as follows: 



1 K^K^M^ , ,2 1 h^K:^M^ , , 



The first term in (2.22) contributes to the 1-loop infiationary potential (2.12), while the 
remaining Q^-independent terms lift the fiatness of the itij-direction. Assuming that <S 
A^ and M '^\S\ towards the end of infiation, we find the well-defined minimum 

{ur) = (2.23) 

We should remark here that the above minimum would remain unaltered, even if the fiat 
direction were a squark or slepton doublet. In this case, only the overall normalization of 



11 



the Q^-independent part of V'^^^u^) would have changed by a factor 1/2. The loop-induced 
VEV of ur generates a squared mass via the Higgs mechanism, which is given by 

This squared mass should be compared with the size of possible negative Hubble- 
induced squared mass terms of order H'^ — k^M^/ (3m|,i), e.g. terms of the form — c| H'^\uji\'^ 
that may occur in V'^^^ (uji) and originate from SUGRA effects. These terms may play some 
role in our model, unless c~if^ < M~^. The latter condition may be translated into the 
inequality 

c- < _J_^^!!^ (2 25) 

As a typical example, let us consider an inflationary scenario, with A = 2k,, k = 10^^ and 
M = 10^^ GeV. In this case, (2.25) implies that < 0.87 h. Hence, although the required 
tuning of the coefficient to fulfill this last inequality may not be significant for the third 
generation squarks and sleptons, it becomes excessive for the first generation, unless a 
minimal Kahler potential is assumed. It should be stressed here, however, that the deepest 
and hence most energetically favoured minimum for all squark and slepton directions is the 
one related to Ir. In other words, given chaotic initial conditions, the fields are most likely 
to settle to minima of quasi-flat directions involving large Yukawa couplings. In this case, 
radiative effects play an important role in the dynamics of flat directions^. 

Let us finally assume that we are in a situation where the Hubble-induced mass terms 
can be neglected, i.e. c^i = as is the case for a minimal Kahler potential, for example. 
Suppose that the loop-induced VEV of the quasi-fiat direction persists throughout the 
coherent oscillatory regime. In this case, the VEV (2.23) gives rise to masses h{uR) — 
nM/ \/6 in the Q-^u-sector, which do not depend on the Yukawa coupling h. Consequently, 
the inflaton- related flelds of mass \/2k M (see Table 3) will have a large decay rate to those 
massive particles, thus creating a non-thermal distribution. This non-thermal distribution 
will in turn induce T-dependent mass terms which can be larger than the expansion rate 
H{T) at some temperature T soon after inflation, such that {ur) will rapidly relax to zero. 
Of course, one might think of contemplating conflgurations where multiple flat directions 
have VEVs which contribute constructively to the masses of both if„ and Hd, such that 
all inflaton and waterfall particle decays would be kinematically forbidden. However, we 



^We should note that the evolution of flat directions during the waterfall and coherent oscillation periods 
is a non-equilibrium dynamics problem. Moreover, no theoretical methods yet exist that would lead to a 
practical solution to this problem, even though effective potential corrections to the flat directions as the 
ones considered here are expected to be relevant during the above cosmological periods. 
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consider such a possibility as a bit contrived. It is therefore reasonable to assume that, 
provided (2.25) is fulfilled, reheating and equilibration of all MSSM degrees of freedom will 
take place in the Fu-teim hybrid model and in all supersymmetric models of inflation that 
include an unsuppressed renormalizable operator of the form SHuHa. 

2.3 Topological Defects and GUT Embeddings 

As we mentioned in the Introduction, topological defects, such as domain walls, cosmic 
strings or monopoles, may be created at the end of inflation, when a symmetry group G, 
local, global or discrete, breaks down into a subgroup H, in a way such that the vacuum 
manifold M — G/H is not trivial. Specifically, the topological properties of the vacuum 
manifold M under its homotopy groups, 7r„(M), determine the nature of the topological 
defects [22,23]. Thus, one generally has the formation of domain walls for 7ro(M) 7^ I, 
cosmic strings for 7ri(M) 7^ I, monopoles if 7r2(M) 7^ I, or textures if 7rn>2(M) 7^ I [22]. 
For example, for the SSB breaking pattern U(l)x — > I in the waterfall sector, the first 
homotopy group of the vacuum manifold is not trivial, i.e. 7ri(U(l)/I) = Z. In this case, 
cosmic strings will be produced at the end of infiation. In general, the non-observation 
of any cosmic string contribution to the power spectrum P]z at the 10% level introduces 
serious constraints on the theoretical parameters of hybrid infiation models. 

A potentially interesting infiationary scenario arises if the waterfall sector possesses 
an SU(2)x gauge symmetry. In this case, the SSB breaking pattern is: SU(2)x — > I, i.e. the 
group SU(2)x breaks completely. It is worth stressing here that this is a unique property of 
the SU(2) group, since the breaking of higher SU(A'^) groups, with N > 2, into the identity I 
is not possible. Moreover, an homotopy group analysis gives that 7ro,i,2(SU(2)x/I) = I, im- 
plying the complete absence of domain walls, cosmic strings and monopoles. The only 
non-trivial homotopy group is 7r3(SU(2)x/I) = Z, thus signifying the formation of tex- 
tures, in case the SU(2)x group is global. If the SU(2)x group is local, however, observable 
textures do not occur. Since their corresponding field configurations never leave the vac- 
uum manifold, the would-be textures can always be compensated by local SU(2)x gauge 
transformations [22] . It is therefore essential that the X-symmetry of the waterfall sector 
is local in the Fo-term hybrid model. 

It is now interesting to explore whether generic scenarios exist, for which the waterfall 
gauge groups U(l)x or SU(2)x of the Fo-term hybrid model may, partially or completely, 
be embedded into a GUT. As a key element for such a model-building, we identify the 
maintenance of iP-parity conservation in the X-gauged waterfall sector, which is discussed 
in detail in Section 4.1. In order to preserve D-parity, the waterfall sector should be 
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somehow 'hidden' from the perspective of the SM gauge group Gsm- This means that the 
SM fields must be neutral under X and vice versa, the X-gauge and waterfall sector fields 
should not be charged under Gsm- Consequently, we have to require, as a GUT breaking 
route, that the waterfall X-gauge group and the GUT-subgroup that contains Gsm factor 
out into a product of two independent groups without overlapping charges. 

It is reasonable to assume that the GUT-subgroup is broken to Gsm before or while 
inflation takes place. Then, possible unwanted topological defects due to the various stages 
of symmetry breaking from the GUT-subgroup down to the SM will be inflated away. 
Notice that we do not have to require that the GUT-subgroup breaking scale is higher 
than the respective X-symmetry breaking scale, but only that the reheat temperature T^^h 
is low enough such that no symmetries of the GUT-subgroup are restored during reheating. 
A related discussion within the context of SO (10) may be found in [36]. 

Let us first investigate whether a 'hidden' gauge group U(l)x related to the waterfall 
sector can be embedded into a GUT. Although 'hidden' U(l)'s naturally arise in models of 
string compactification [37], our interest here is to identify possible U(l)x factors that can 
be embedded into a simple GUT. Given the above criterion, the frequently discussed GUT 
based on SO(IO) should be excluded, since it does not contain 'hidden' U(l)x groups [38]. 
As a next candidate theory, we may consider the exceptional group E(6), with the SSB 
breaking path E(6) U(l) x SO(IO). The fundamental representation of E(6) is the chiral 
27 p representation, which branches under U(l) x SO (10) as follows: 

27 F = (4, 1) + (-2, 10) + (1, 16) . (2.26) 

Although the SM particles may fit into 16, they are not neutral under the extra U(l). 
Higher representations, such as (0,45) stemming from 78 of E(6), are neutral under the 
U(l) factor, but they are not suitable to properly accommodate all the SM particles. 

We therefore turn our attention to possible breaking patterns of maximal groups that 
contain a 'hidden' SU(2)x factor. A promising example is E(6) D SU(2)x X SU(6), where 
the fundamental representation 27 p follows the branching: 

27 F = (2,6) + (1,15) . (2.27) 

Under SU(6) D SU(5) x U(l), 15 is an antisymmetric representation of SU(6) and one of 
its branching rules is 

15 = (5,-4) + (10,2) . (2.28) 

However, we need a 5 of SU(5), together with 10 in (2.28), in order to appropriately describe 
all SM fermions. This shortcoming may be circumvented by adding an extra 27 p of E(6) to 
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the spectrum, where the missing 5 may be obtained from the complex conjugate branching 
of (2.28). Such an extension of the particle spectrum may even be welcome to resolve the 
proton stability problem, through a kind of split multiplet mechanism [39]. Within the 
framework of SUSY, the quark and lepton Yukawa interactions may be generated via the 
introduction of a pair of the multiplets 27^, 27 h- Finally, in such an E(6) unified scenario, 
the 3 right-handed neutrinos can only appear as singlets. 

Another possible GUT scenario that complies with our criterion of a hidden SU(2)x 
is E(7) D SU(2)x X S0(12). The fundamental representation is 56^ and branches under 
SU(2)x X S0(12) as follows: 

56j7 = (2,T2) + (1,32). (2.29) 

Subsequently, S0(12) breaks spontaneously into SO(IO) x U(l), where 32 = (16,1) + 
(16, —1) is a vector-like representation. However, one may well envisage a string-theoretic 
framework, in which orbifold compactification projects out the undesirable anti-chiral 
states. Then, all SM particles, including right-handed neutrinos, will be contained in one 
of the 16's of 32. Related discussion of missing or incomplete multiplets due to orbifold 
compactification may be found in [40]. 

Building a realistic GUT model from the blocks stated above lies beyond the scope 
of this paper. We have demonstrated here, however, that the embedding of an SU(2)x 
gauge group into a GUT, which is hidden but nevertheless takes part in the gauge coupling 
unification, appears feasible within E(6) and E(7) unified theories. 

We conclude this section by observing that the presence of the singlet inflaton field S 
offers alternative options, for suppressing the heavy Majorana neutrino masses within SUSY 
GUTs. As an example, we mention the breaking scenario, where SO (10) — > SU(5) via the 
VEV of a 126i^ Higgs representation and the usual superpotential term 16^(126^^) 16^ 
induces heavy Majorana masses of the GUT scale Mqut- Given that the above renormal- 
izable operator is forbidden by some i?-symmetry, the presence of an i?-charged inflaton S 
may give rise to a drastic suppression of the GUT-scale Majorana mass, through a superpo- 
tential term of the form S 16i7'(126/f) 16F/mpi. Since S receives a VEV of order Msusy//* 
in general F-term hybrid models [cf. (2.4)], one naturally obtains heavy Majorana neutrino 
masses of order Msusyj if ~ {12Qh) /mpi ~ 10~^. Such values of k do satisfy the current 
inflationary constraints which we discuss in the next section. 
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3 Inflation 



Here, we first briefiy review in Section 3.1 the basic formafism of inflation, including the 
constraints from the non-observation of cosmic strings in the power spectrum P-ji of the 
CMB data. Then, in Section 3.2, we present our numerical results for two scenarios: (i) the 
minimal SUGRA (mSUGRA) scenario and (ii) the next-to-minimal SUGRA (nmSUGRA) 
scenario. In particular, we exhibit numerical predictions for the spectral index Ug and 
discuss its possible reduction in the nmSUGRA scenario. Finally, we analyze the combined 
constraints on the fundamental theoretical parameters A, p, and M, which result from 
the recent CMB observations and inflation. 



3.1 Basic Formalism 

According to the inflationary paradigm [6] , the horizon and flatness problems of the stan- 
dard Big-Bang Cosmology can be technically addressed, if our observable Universe has 
undergone an accelerated expansion of a number 50-60 of e-folds during inflation. In the 
slow-roll approximation, the number of e-folds. He, is related to the inflationary potential 
through: 

1 /•'Aexit y. 

-r #7^ - 55. (3.1) 

Hereafter, a prime on Vinf will denote differentiation with respect to the inflaton field 
(f) = y^KeS. In addition, 0exit is the value of 0, when our present horizon scale crossed 
outside inflation's horizon and 0end is the value of (p at the end of infiation. In the slow-roll 
approximation, the field value 0end is determined from the condition: 

max{e(0end), |??(0end)|} = 1, (3.2) 

where 

We have checked that the slow-roll condition (3.2) is well satisfied up to the critical point 
0end = V^M, beyond which the waterfall mechanism takes place. We also find that the 
slow-roll condition remains valid, even within the nmSUGRA scenario with Ch ^ and 
with appreciable non-renormalizable SUGRA effects. Finally, we note that the assumed 
value of A/"e — 55 in (3.1) is slightly higher than the one computed consistently from (5.17), 
which is about 50 for our low-rchcat cosmological scenario. However, our numerical results 
concerning Pji and do not depend on such a 10% variation of J\fe in any essential way. 
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The power spectrum P-ji is a cosmological observable of the curvature perturbations, 
which sensitively depends on the theoretical parameters of the inflationary potential. The 

1 /2 

square root of the power spectrum, , may be conveniently written down as 

pl/2 ^ 1 ^nf^(0exit) a\ 

^ 2^7rm|i |V^:f(</'exit)| ■ ^' ' 

The recent WMAP [2, 4] results, which are compatible with the ones suggested for the 
COBE normalization [1], require that 

P^/^ ~ 4.86 X 10"^ . (3.5) 

In addition to scalar curvature perturbations, tensor gravity waves and cosmic string 
effects may also contribute to P-jz. In the F^-term hybrid model with an Abelian U{l)x 
waterfall sector, cosmic strings arise after the SSB of the gauge symmetry (see also our 
discussion in Section 2.3). In this case, additional constraints are obtained from the 
non-observation of cosmic string effects on Pji [41,42]. The evaluation of such effects in- 
volves a certain degree of uncertainty in the numerical simulations of string networks [43]. 
Nevertheless, the common approach taken to cosmic string effects [29, 44] is to require 
that their contribution {Ptz)cs to the power spectrum P-ji does not exceed the 10% level, 
i.e. {Pn) cs / Pn ~ 0.1. In detail, we require that 

(p4^')cs < 1.54xl0-^ (3.6) 

The cosmic string contribution {Pn)cs to the power spectrum may be computed by 



/.es = 2nMhM, UP) ^ \ . . . 2 (3-8) 



(PnX - yes , (3.7) 

where the tension of the cosmic strings, Hcs, is calculated using the formulae: 

1.04 , for /? > 10^2 

2.4/ln(2//?), for (3 < 10 

In (3.8) the argument /3 is given hy P — k^/ {'2g^), while the U(l)x gauge coupling constant 
g is considered to assume the value g ~ 0.7 as is the case in GUT models. The central value 
of the parameter yes is 8.9 and its error margin lies in the interval [6.7,11.6], according to 
the analysis in [42]. 

The recently announced three-years results of WMAP [4] improved upon the precision 
of a number of other cosmological observables. The merits of an inflationary model can be 
judged by comparing its predictions for the scalar spectral index, rig, the tensor to scalar 
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ratio, r, and the running of rig, dUs/dlnK,, with the CMB data. In the Fo-term hybrid 
model, r = 16e((6cxit) is much lower than the WMAP bound, i.e. well below 10~^, and 
dUg/dlnK is always smaller than 10^^ and so unobservable. In addition, the spectral index 
Ug in our model may well be approximated as follows: [6] 

Us = l-6e(0exit) + 277(0exit) 1 + 277(0exit), (3.9) 

since e is neghgible. The predicted value needs to be compared with the recent WMAP 
results [4]: 

= 0.9511^^9 ■ (3.10) 
The latter is translated into the double inequality, 

0-913 < < 0.981 , (3.11) 

at the 95% confidence level (CL). 

The result (3.11) brings under considerable stress minimal F-term hybrid inflation 
models [9]. This is due to the fact that these models predict Ug extremely close to unity 
without much running. To be precise, when the radiative corrections dominate the slope 
of the potential, we obtain 

ris ~ 1 - 1/K - 0.98, (3.12) 

for A/'e = 55. On the other hand, if the non-rcnormalizablc operator \S\'^ in Vsugra of (2.13) 
dominates the slope of the potential of a mSUGRA model with ch = [12], we obtain a 
blue-tilted spectrum, with 

' + ^1, - 2APK ' ' ■ 

A possible Hubble-induced positive term +c^iJ^|S'p in Vsugra [27,29] implies an even 
more pronounced blue spectrum and is therefore excluded by the current WMAP data. 

As noticed earlier in [45] and elaborated further in Ref. [46], agreement of theory's 
prediction for ng with observation strongly suggests the presence of a negative Hubble- 
induced mass term —c'jjH'^\S\'^ in Vsugra, thereby clearly disfavouring the minimal Kahler 
potential. In our analysis, we therefore consider the following next-to-minimal form for the 
Kahler manifold [27]: 

Ks ^ \S\' + ksP^ , (3.14) 

where the constant ks can be either positive or negative. Substituting (3.14) into the 
general formula for the F-term type contributions to the SUGRA potential (see, e.g. [8]), 



nKs')lFj - 3 



\W 



2 1 



mil 



(3.15) 
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we arrive at the result (2.13) with = 3ks, after neglecting higher-order terms that are 
small for \ch\ ^ 0.2. In (3.15), are the SUGRA-generalized F-terms and {Kg^)l is the 
so-called inverse metric of the Kahler manifold, where a superscript (subscript) index i or 
j on Ks denotes differentiation with respect to S (S*). 

The aforementioned nmSUGRA inflationary potential, with a negative Hubble-induced 
mass term, reaches a local minimum and maximum at the points (pmin and 0max; respec- 
tively. These points can be estimated by 



(gkW + 12A2 + 9p2)'^', 0n.in ^ y|cHmpi . (3.16) 



Attch 

For relevant parameter values, for which ^max < <Pm.m, and under convenient initial condi- 
tions, the so-called hilltop inflation [45] can take place, where rolls from 0max down to 
smaller values, such that 0exit < 0max- In this nmSUGRA scenario, the value of Us can be 
significantly lowered and can be approximately given by 

ris ~ 1 - - 4 . (3.17) 

As we will show more explicitly in the next section, the spectral index Ug can be easily 
driven into the range of (3.11), for values of ch ~ 0.1. The presence of the second next-to- 
minimal term proportional to ks in (3.14) modifies the analytic expressions of (3.1), (3.3) 
and (3.5) [46]. However, these modifications turn out to be numerically insignificant for 
the predicted values of Afg and P-ji, if ch is not very large, e.g. ch ~ 0.2. 



3.2 Numerical results 

In our numerical estimates, we use the full expression for the inflationary potential Vinf 
given in (2.8), which consists of the tree-level, 1-loop and SUGRA contributions, given 
in (2.9), (2.12) and (2.13), respectively. We will ignore all soft SUSY-breaking terms, but 
the tadpole term as- To facilitate our numerical analysis, we introduce the real tadpole 
parameter a.s, which is defined, in terms of the Lagrangian parameter as, by the relation: 

a,s — — 2|a5| cos (argos -I- argS") . (3.18) 

For any given value of k, A, p, a.s and ch, we determine 0exit and M, by imposing the 
conditions (3.1) and (3.5) for the number j\4 of e-folds and the power spectrum P^'^, 
respectively. In addition, we compute Ug by means of (3.9). Our results are presented in 
Fig. 1 for the mSUGRA scenario and in Fig. 3 for the nmSUGRA scenario. They will be 
analyzed in more detail in the following two subsections. 
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Figure 1 : The values of the inflationary scale M allowed by (3.1 ) and (3.5) (a) and the predicted 
values of the spectral index Ug (b) as a function of k for J\f = 1 and p = X = k (light grey lines) 
or p = \ = Ak (grey lines), including the one-loop radiative corrections (dashed lines) or the 
mSUGRA (ch = 0) contributions with as = I TeV (solid lines). The upper bound of (3.6) for 
Ucs = 6.7, 8.9, 11.6 (from top to bottom) [cf (3.10)] is also shown by thin lines (a) [(b)]. 
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Figure 2 : The allowed values ofX/K versus p/n for the mSUGRA scenario with M = 2 x 10-^^ GeV 
and K = 0.05 (dark grey line), k = 0.01 (grey line) or k = 0.005 (light grey line). 

3.2.1 The minimal SUGRA scencirio 

Here, we present numerical results for the mSUGRA scenario. The values of the inflationary 
scale M allowed by (3.1) and (3.5) and the predicted values of rig, as functions of k, for 
p — X — K, (light grey hues) and p — X — Ak, (grey lines), are displayed in Fig. 1(a) and 1(b), 
respectively. Dashed lines indicate results obtained, when only the 1-loop contribution 
to Vinf is considered and as is set to zero, whilst solid lines represent numerical values 
obtained, if the remaining contributions arc included, namely those coming from (2.9) with 
a5 = 1 TeV and (2.13) with Ch = 0. In Fig. 1, we observe that as the common value 
for p, X and k, increases, M and Us increase as well. In particular, M gets closer to the 
GUT-scale value 2 x 10^^ GeV for k ~ 10"^, unhke the case A = p = 0, where M takes on 
much smaller values at this point [9,12,29]. 

It is now not difficult to identify in Fig. 1 the regimes, in which the different contri- 
butions to Vini dominate. More explicitly, for k > 4 x 10~^ and p — X — k,oj:k> 10'^ and 
p — X — Ak, the non-renormalizable SUGRA term of (2.13) dominates and drives Us to val- 
ues close to or larger than 1 [cf. Fig. 1(b)]. On the other hand, for 4 x 10~^ ^ ^ 4 x 10~^ 
and p = A = «;or4x 10~^ ^ i$ 10~^ and p = X = 4k, the 1-loop corrections (2.12) 
dominate, in which case the spectral index rig takes on the predicted value ~ 0.98 given 
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M = 2 X 10^^ 
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K - 


= 0.005, M = 


10^6 GcV 


K 


0cxit 


Us 


Ch 


0min 


0max 


0cxit 




0.005 


6.28 


1.017 


0.07 


8.4 




5.1 


0.978 


0.001 


2.14 


0.99 


0.14 


16.5 


10.5 


8.1 


0.955 


0.0005 


1.51 


0.99 


0.18 


21.7 


12.8 


11.2 


0.941 



Table 1 : The values of <pcxit (in units \piM) and Ug for several k 's along the curves in Fig. 2 and 
the values of ^min; ^max, ^^'exit (^^ units y/2M) and Us for several ch's along the curves in Fig. 5. 

in (3.12). Finally, for /t < 4 x 10~^ and p = \ = Koip = \ = 4k,, the tadpole term in (2.9) 
starts playing an important role. As M increases, the non-renormalizable SUGRA term 

1 /2 

of (2.13) becomes again important [12,29]. In this case, the prediction for and is 
almost independent of p and A, as expected. For lower values of ag, the solid lines in the 
latter regime would eventually approach the dashed hues [29]. In Fig. 1(b), we also indicate 
with a thin line the 95% CL upper hmit on Us stated in (3.11). Clearly, a mSUGRA version 
of the Fp-term hybrid model appears to be disfavoured by the most recent WMAP results. 

In addition, we show in Fig. 1(a) upper limits due to cosmic string effects based 
on (3.6), for ?/cs = 6.7, 8.9, 11.6 (from top to bottom). Such constraints are only relevant 
for an Abelian realization of the waterfall-gauge sector. We observe that the presence of 
cosmic strings severely restrict the available parameter space of the U(l)x F^-term hybrid 
model. As we discussed in Section 2.3, however, these constraints do no longer apply, if 
the waterfall-gauge sector realizes an SU(2)x gauge symmetry. Since the dimensionality 
of the representation is jV' = 2 in this case, the allowed range of M as a function of k, 
slightly changes. In fact, the allowed values of M become marginally larger than the ones 
already shown in Fig. 1(a) by up to 12%, ior p — X — while they stay at the 1% level, for 
p — X — Ak. Likewise, the predicted values of Ug remain almost unaffected at the 2% level, 
from those presented in Fig. 1(b). Obviously, as p and A gets larger than the difference 
between the M —1 and M — 2 cases becomes practically unobservable. 

Finally, in Fig. 2 we plot the allowed values oi X/k versus p/n, subject to the con- 
straints (3.1) and (3.5), for M = 2 x 10^^ GeV (close to the GUT scale) and for different 
values of k — 0.005 (dark grey line), k — 0.001 (grey line) k — 0.0005 (light grey line). 
Along these contour lines, 0exit and Us remain constant and equal to the values presented 
in Table 1. We observe that as k increases, 0exit and Us increase as well. 



22 





K 

Figure 3: The values of the inhationary scale M allowed by (3.1 ) and (3.5) (a) and the predicted 
values of the spectral index Us (b) as a function of k for J\f = \ and p = X = k (light grey lines) 
or p = \ = Ak (grey lines), for the nmSUGRA scenario with cjj = 0.07 (dashed lines) or 
ch = 0.14 (solid lines). In both cases we take as = 1 TeV. The upper bound given in (3.6) (for 
yes = 6.7, 8.9, 11.6 from top to bottom) [cf. (3.10)] is also depicted by thin lines (a) [(b)j. 
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3.2.2 The next-to-minimal SUGRA scenario 

We now turn our attention to the nmSUGRA scenario. Ahhough we take the tadpole term 
to be a5 = 1 TeV, its impact on our results turns out to be insignificant for the whole range 
of parameters we have scanned. The values of the infiationary scale M allowed by (3.1) and 
(3.5) and the predicted rig as a function of k are presented in Fig. 3(a) and 3(b), respectively, 
ior p — X — K, (fight grey fines) and p — X — Ak, (grey fines). We consider the two cases: 
ch — 0.07 (dashed lines) and ch — 0.14 (sofid fines). As in the case of mSUGRA, M and 
Us increase, with increasing p, X and k. Moreover, as k decreases, the non-renormalizable 
SUGRA contribution in (2.13) becomes subdominant and Us decreases. Such a reduction 
becomes even more drastic with increasing ch, as can be easily inferred from Fig. 3(a), 
where the 95% CL upper bound on Us [cf. (3.11)] is indicated with a thin horizontal line 
on the same plot. In stark contrast to the mSUGRA scenario, we observe that our model 
can become perfectly consistent with the recent WMAP result for 0.04 < 0.22. Note 
that the various lines terminate at large values of for which the two restrictions (3.1) 
and (3.5) cannot be simultaneously met. 

It is interesting to further investigate the inflationary dynamics described by Vinf in 
the presence of a negative Hubble-induced mass term. To this end, we exhibit in Table 2 
the values of Ch-, 0mm; 0max! 0cxit (i^ uuits ^/2M) and the inflationary scale M (in units 
of 10^^ GeV) which are obtained for different values of k, assuming that X = p = k 
OT X = p = 4k, and for flxed values of rig, i.e. rig = 0.913, 0.951, 0.981, compatible 
with the 95% CL limits given in (3.11). In addition, we present values for the parameter 
Aexit = {4>ma.x — 0cxit ) / 0cxit , which somehow quantifies the degree of tuning required in the 
initial conditions of inflation. The entries without a value assigned (in Tables 1 and 2) 
mean that the respective inflationary potential Vinf has no distinguishable nearby local 
maximum 0max- We notice from Table 2, that as Ug decreases with flxed values of k,, Ch 
increases while M and Agxit decrease. Moreover, for flxed values of and decreasing n, 
Ch and M decrease and 0exit approaches 0max- On the contrary, with increasing k, X and 
p, the inflationary scale M increases and the parameter Agxit becomes larger. We have 
checked that the inequality 0max > <Pexit is fulflUed along the lines presented in Fig. 3. In 
this respect, we also note that 0min is in general much larger than 0max especially for low 
values of Ug. 

It is important to observe from Table 2 that there is a degree of tuning required for 
the values 0cxit with respect to 0max- For values of k > 10^^, we find that the degree of 
tuning required is not very serious, i.e. Acxit ~ 10%. However, the situation becomes rather 
delicate as k, gets smaller than 10~^, for Ug ^ 0.97. In this case, we flnd that 0max ~ 4>exit, 
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19.0 
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0.058 


1.30 


4.7 
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0.0005 


0.178 


0.46 


54 


1.53 


1.49 


0.028 


0.129 


0.68 


26 


1.45 


1.33 


0.08 


0.070 


1.00 


9.6 


1.83 


1.30 


0.29 



Table 2: The values of ch, M (in units 10^^ GeVj (f)min, 0max, ^z!>exit (in units \piM) and Aexit = (0max - ^z!'exit)/^max, for selected 
values of k, A and p, and for fixed values of the spectral index rig. 




Figure 4: The parameter values {k,ch) allowed by (3.1), (3.5) and (3.10) in the nmSUGRA 
scenario, for p = X = k (light grey hatched area) and p = X = An (grey hatched area). The grey 
(light grey) line has been obtained by fixing Ug to its central value given in (3.10), for p = X = An 
(p = X = k). 



leading to a substantial tuning at the few per cent level in the initial conditions of inflation. 

As in the mSUGRA case, we also show in Fig. 3(a) the upper bounds resulting 
from cosmic-string effects [cf. (3.6)], for t/cs = 6.7, 8.9, 11.6 (from top to bottom). As 
mentioned above, these constraints are only relevant for an Abelian waterfall-gauge sector 
with dimensionality M — 1. However, unlike in the mSUGRA case, these restrictions 
appear less harmful, since the inflationary scale M assumes smaller values (see also Table 2) 
and the tadpole term becomes unimportant. Thus, larger values of k up to order 10~^ can be 
tolerated in this case. For the non-Abehan SU(2)x Fp-term hybrid model, the restrictions 
from considerations of cosmic-string effects are totally hfted and the lines depicted in Fig. 3 
only vary within the few per cent level. Such a variation becomes even smaller li p > k 
and/or X> k. 

In Fig. 4, we present the parameter space {k,,ch) which is allowed by the condi- 
tions (3.1), (3.5) and (3.10) in the nmSUGRA scenario. The light grey (grey) hatched area 
indicates the allowed region ior p — X — Ak (p — X — k). The lower (upper) boundaries of 
the allowed regions correspond to the upper (lower) bound on Ug, cf. (3.11), while the solid 
hues correspond to the central value of Ug, cf. (3.10). We find that values of ch ~ 0.2 and 
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Figure 5: The allowed values of p/k versus \/k for the nmSUGRA scenario with k = 0.005, M = 
10^^ GeV and ch = 0.18 (dark grey line), ch = 0.14 (grey line) or ch = 0.07 (light grey line). 

K ~ 0.05 are still possible in a nmSUGRA extension of the Fo-term hybrid model. 

Finally, we plot in Fig. 5 the allowed values of A//t versus p/n, on account of the 
inflationary constraints (3.1) and (3.5), for k = 0.005, M = 10^^ GeV, and for = 0.18 
(dark grey line), Ch = 0.14 (grey line) and Ch = 0.07 (light grey line). We have selected 
a slightly lower value for M, because no viable nmSUGRA scenarios seem to exist with 
acceptable values for rig, if M = 2 x 10^^ GeV and ch > 0.07. Along the contour lines 
in Fig. 5, 0mm, 0max, 4>exit ^ud Ug remain constant and equal to the values presented in 
Table 1. We observe that as ch increases, 0exit approaches 0max, 0min increases, while rzg 
decreases. This kinematic behaviour is in agreement with our discussion related to (3.16) 
and (3.17). 



4 Preheating 

As stated in the Introduction, gravitinos, if thermally produced during the early stages 
of the evolution of the Universe, will spoil the successful predictions of BBN [47]. Their 
disastrous consequences may be avoided, if the reheat temperature T^eh of the Universe is 
not very high. In fact, depending on the decay properties of the gravitino, it should be 
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Trch ~ 10^-10^° GeV [51,56]. This fact leads to a tension between the allowed range of Trch 
and the natural scale of hybrid inflation Af , which is of order ~ 10^^ GeV. The traditional 
way taken to get around this problem is to consider scenarios where the decay rate of the 
inflaton to SM particles is extremely suppressed, e.g. by suppressing all possible couplings 
of the inflaton to the SM flelds. 

In this and next sections, we present in detail an alternative solution to the above 
gravitino overabundance problem [13]. Our solution relies on the huge entropy release 
caused from the late out-of-equilibrium decays of the supermassive waterfall particles. The 
entropy produced through this mechanism is sufficient to reduce the gravitino abundance 
Yq to levels compatible with BBN limits discussed in detail in Section 5. Figure 6 gives 
a schematic representation of the post-inflationary dynamics of the early Universe, as is 
predicted by the Fo-term hybrid model. Shortly after inflation ends, the energy density 
of the Universe is predominantly stored to coherently oscillating inflaton condensates 
which scale as a~^, where a is the usual cosmological scale factor describing the expansion of 
the Universe. The coherent oscillations of the inflaton-related condensates also give rise to 
another non-perturbative mechanism called preheating. During preheating, waterfall gauge 
particles of energy density Pg are produced almost instantaneously, which are absolutely 
stable if a £)-parity, an analogue of the usual i?-parity in the MSSM, is conserved. Then, 
the following scenario visualized in Fig. 6 emerges. First, Pg/p^ remains constant during the 
epoch of coherent oscillations, since both pg and p^ behave as matter energy densities and 
scale as during this period. The constancy of pg/ p^ ceases to hold, when the coherently 
oscillating inflaton condensates decay and their energy density p^ gets distributed among 
light relativistic degrees of freedom. As a consequence of the latter, p,^ will be oc a~^, 
whilst Pg will still be oc a~^. If the initial value of Pg/pn is not very suppressed, e.g. it is of 
order 10~*^-10~^, the waterfall gauge particles will eventually dominate the energy density 
of the Universe, leading to a second matter dominated epoch which will last until these 
particles decay via D-parity violating couplings. This is expected to produce an enormous 
entropy release and so reduce the gravitino-to-entropy ratio to values compatible with 
BBN constraints. 

The discussion in this section is organized as follows: in Section 4.1, we pay special 
attention to D-parity and derive the particle spectrum of the combined inflaton-waterfall 
sector in the supersymmetric limit of the theory. In addition, we compute the decay rates of 
all inflaton-related and waterfall gauge particles. Finally, in Section 4.2, we discuss how the 
waterfall gauge particles are instantaneously produced through preheating and calculate 
the resulting energy density pg carried by these particles. 
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Figure 6: Schematic representation of the thermal history of the Universe in the Fu-term hybrid 
model. 

4.1 D-Parities and the Inflaton- Waterfall Sector 

Let us first consider a model with a U(l)x gauge-symmetric waterfall sector. The case of 
a waterfall sector realizing a non-Abelian SU(2)x gauge symmetry is analogous and will 
be discussed later. In terms of superfields, the minimal gauge-kinetic Lagrangian of the 
U(l)x model reads: 

where Vx is the U(l)x vector superfield and Wa {Wo) are their respective chiral (anti- 
chiral) field strengths. The latter are given by 

Wa = - — D^e-^s^^Dae^'^^), = — D\e^<^^^ D^e-^^^^) , (4.2) 

where and are the usual SUSY-covariant derivatives which are irrelevant for our 
discussion here. The minimal gauge-kinetic Lagrangian (4.1) possesses the discrete sym- 
metry 

D: X, ^ X2, Vx 'Vx, (4.3) 
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whereas all other superfields do not transform. It is not difficult to verify that the complete 
Fo-term hybrid model, including the superpotential (2.1) and its associated soft SUSY- 
breaking sector, is invariant under the discrete symmetry (4.3) in the unbroken phase of the 
theory. After the SSB of U(l)x, the waterfall fields receive the VEVs: (Xi) = (X2) = M. 
Thus, the above discrete symmetry survives even in the spontaneously broken phase of the 
theory. Since the discrete symmetry acts on a gauged waterfall sector, it manifests itself as 
a kind of parity, which we call i?-parity. 

It therefore proves convenient to choose a weak basis where the fields are eigenstates of 
D-parity. To this end, we define the linear combinations in terms of the waterfall superfields 

X± = -^(Xi ± X,) . (4.4) 

Evidently, the superfield X+ i^-) has even (odd) D-panty; its D-parity quantum number 
is +1 (—1). The vector superfield Vx, which is already a £)-parity eigenstate, has odd 
D-parity. All remaining fields, including the infiaton superfield S and the other MSSM 
superfields, have positive D-parity. 

As a consequence of D-parity conservation, all D-odd particles will be stable, in as 
much the same way as the usual i?-parity guarantees that the LSP of the MSSM is stable. 
As we explicitly mentioned in Section 2.1, the simplest way to break £)-parity is to add a 
FI D-term to the model, e.g. 

Ui = -l^li j d'OVx = -^ml,D, (4.5) 

where D is the auxiliary component of the vector superfield Vx- It is obvious that Cpi flips 
sign under the discrete symmetry (4.3). Other mechanisms of explicitly breaking D-parity 
are discussed in Appendix A. 

We now calculate the particle spectrum of the infiaton- waterfall sector in the presence 
of a subdominant FI D-term mpi and in the supersymmetric limit of the theory. With this 
aim, we expand the scalar D-parity eigenstates X± about their VEVs: 

X± = (X±) + 1=(^R^ + i/±) . (4.6) 

The VEVs (V±) are determined from the minimization conditions of the combined F- and 
D-term scalar potential 

Vfd = F^Fs + \d\ (4.7) 
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where 

Since SUSY is preserved after the SSB of U(l)x, the scalar potential Vpo will vanish at its 
ground state, i.e. (Vpo) = 0. Consequently, to leading order in mpi/M, the VEVs of the 
scalar inflaton-waterfall fields are 

(5) = 0, = V2M, {X_) = , (4.9) 

where v — m|j/ (2M) . Notice that the VEVs of the F- and D-terms vanish through order 
mpi/M considered, i.e. (D) = and (F5) = 0{m^jM^). 

To derive the mass spectrum, we expand the potential about its ground state up 
to terms quadratic in all the fields involved. We first consider the F-terms. To order 
v/M {— m|j/M^), we find the approximate mass eigenstates: 

S^^fcp + ia), R+--^R-, (4.10) 

All the above fields, consisting of 4 bosonic degrees of freedom in total, share the common 
mass 

m„ = V2kM . (4.11) 

As a consequence of SUSY, the corresponding 4 fermionic degrees of freedom form a Dirac 
spinor ip^, which also has the same mass (4.11). We refer to these particles as inflaton- 
related or re-sector particles. 

The remaining scalar fields receive their masses from the D-term of the scalar potential 
VpD in (4.7). Performing an analogous calculation as outlined above, we obtain to order 
v/M the scalar mass eigenstates: 

/_ + —/+, R_ + —R: . (4.12) 
2M +' 2M + ^ ^ 

The first field is absorbed by the longitudinal component of the U(l)x gauge field V^, 
via the Higgs mechanism. In the supersymmetric limit, all these fields, which mediate 4 
bosonic degrees of freedom, are degenerate and characterized by the common mass 

rUg = gM . (4.13) 

Like in the re-sector case, the respective 4 fermionic degrees of freedom will make up a 4- 
component Dirac spinor of mass rUg. We refer to this group of particles as waterfall gauge 
or gi-sector particles. In Table 3, we present a summary of all the inflaton-related (re-sector) 
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Sector 


Boson 


Fermion 


Mass 


Inflaton 
(/t-scctor) 
/^-parity: +1 


+ 2M ' 
T '" T 




I J 


1 


, /Ou- A/T 


U(l)x 

Waterfall Gauge 
((yf-sector) 
D-parity: —1 


^- + 




y -iA' J 




gM 



Table 3: Particle spectrum of the inSaton and the U(l)j^ waterfall-gauge sectors after inflation, 
where the approximate D -parity for each sector is displayed. The field denotes the U(l) gauge 
boson and A its associate gaugino. The would-be Goldstone boson related to the longitudinal 
degree of appears in the square brackets. 



and waterfall-gauge (g'-sector) particles. As can also been seen from the same Table 3, re- 
sector particles are predominantly D-even, whereas the g^-sector ones have approximately 
D-odd parity. 

It is now interesting to calculate the decay rates of the k- and gf-sector particles and 

analyze their implications for the reheat temperature of the Universe. Starting with the 

singlet field S, it decays predominantly into pairs of charged and neutral higgsinos, h^di 
~o 

^u,d-> ^u,di i^^o pairs of right-handed Majorana neutrinos 1^1,2,3 r- On the other hand, 
the scalars R+ and /+ decay into the SUSY-conjugate partners of the aforementioned fields 
at the same rate. In fact, we find a common decay rate for each of the re-sector particles: 

r. = ^ + 3p^)m«. (4.14) 

The reheat temperature resulting from these perturbative decays of the re-sector particles 
may be estimated using the relation = H{Ti^), where the Hubble parameter H{T) is 
given in the radiation dominated era of the Universe. In this way, we obtain 

where — 240 is the number of the relativistic degrees of freedom in the Fp-term hybrid 
model. Substituting (4.14) and (4.11) into (4.15), we arrive at the expression: 

, / , . N 1/2 

7; = 8.1 • 10^^ GeV X 



re(4A2 + 3p2) 



inieOeV 



(4.16) 
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Assuming that no relevant amount of entropy is released during the subsequent thermal 
history of the Universe, the gravitino constraint on the reheat temperature ^ 10^ GeV 
requires that each individual coupling k., A and p must be smaller than about 10^^, if 
M ~ 10^^ GeV. Further details are given in Section 5. 

The above unnatural tuning of all inflaton couplings to SM fields may be avoided, if 
the large entropy release from the late decays of the (^-sector particles is properly considered. 
An extensive discussion of this issue is given in Section 5. Here, we simply compute the 
decay rates of the (jf-sector particles which are induced by a non- vanishing FI term mpi. 
In fact, the relevant interaction Lagrangian is given by 

Ant = '-^R-iRl + n)- (4-17) 

As mentioned above, this induces a decay width for the D-odd particle which is easily 
calculated to be 

In close analogy with the K-sector, each g^-sector particle decay rate is equal to F^. 

Let us now consider a model with a waterfall sector based on the SU(2)x gauge 
group. As was mentioned in Section 2.1, the waterfall superfields Xi and X2 are chosen to 
belong in this case to the 2-component fundamental and anti-fundamental representations 
of SU(2)x, respectively. Although the two representations are equivalent for the SU(2) 
case, we nevertheless use this convention, such that its generalization to SU(A^) theories, 
with N > 2, is straightforward. The superpotential is almost identical to the one given 
in (2.1), with the obvious substitution: X1X2 X'[X2. Extending (4.1) to the SU(2)x 
case, the minimal gauge-kinetic Lagrangian is written down 

^Tr {W^W^)6^'\e) + iTr(W^W")5(2)(^) 



-^kin — J ^ 

+ xle'''^^Xi + Xle-'^'^^X 



(4.19) 



In the above, Vx = T" is the SU(2)x vector superfield and W„ = {W^ = T'^) 

are the corresponding non-Abehan chiral (anti-chiral) field strengths in the so-called Wess- 
Zumino (WZ) gauge. The superscript 'T" on Vx, i.e. Vx, indicates transposition that acts 
on the generators = |t" of the SU(2) group, where r^'^'^ are the usual Paufi matrices. 
Finally, the trace in (4.19) is understood to be taken over the group space. 

The minimal SU(2)x gauge-kinetic Lagrangian is invariant under the discrete trans- 
formations 

D,: X, ^ X2, Vx ^ -V^ . (4.20) 
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Notice that under the action of Di in (4.20), the field strengths transform as: Wa — 
— {Wa)'^ and Wa —iWo)'^ in any SUSY gauge, including the WZ gauge. If all other 
superfields do not transform, the complete Lagrangian of the non-Abelian Fo-term hybrid 
model will be invariant under the discrete transformation (4.20) in the unbroken phase of 
the theory. 

Our discussion so far has made no reference to the specific properties of SU(2)x and 
so applies equally well to any SU(A^ > 2) theory. However, in the SU(2)x case, the Fci-term 
hybrid model exhibits an additional Abelian or diagonal discrete symmetry. This may be 
defined by 

L>2: ^ T^Xi, ^2 ^ T^X2, Vx ^ t^Vxt\ (4.21) 

whereas all other superfields do not transform. It is then easy to see that (4.21) implies: 
Wa — ^ Wa in any SUSY gauge and likewise for W a- Since r'^ = {t^Y ^'^d (r'^)^ = I2, 
the invariance of the Lagrangian (4.19) and of the whole model under the action of D2 is 
evident."^ 

We now proceed to compute the particle spectrum of the non-Abelian F^-term hybrid 
model after the SSB of SU(2)x. For this purpose, it is useful to introduce the notation 




(4.22) 



where Z is & generic SU(2)x-doublet or anti-doublet (conjugate) field. The left super- 
scripts ± on Z denote the eigenvalues of the discrete symmetry transformation operator 
D2 = defined in (4.21), and they should not be confused with the corresponding eigen- 
values of the isospin operator of the SU(2)x group. In the unitary gauge, the minimum 
of the scalar potential occurs for the field values 

+Xi = +X2 = M, -Xi = -^2 = 0. (4.23) 

Consequently, the discrete symmetries Di and D2 given in (4.20) and (4.21) remain intact 
after the SSB of the SU(2)x gauge group. Since they act on a gauged waterfall sector, they 
are actually parities. We refer to them as Di- and D2-parities, or collectively as iP-parities. 

^In general, for a waterfall-gauge sector based on an SU(A'' > 2) group, there are N distinct discrete 
symmetries. The first is given by (4.20), while the remaining N — 1 symmetries result from replacing 
with Dn = diag (1, 1,. ..,1,-1,1,. ..,1). The entry —1 occurs at the n position of the A^-dimensional 
diagonal matrix D„, with the restriction 1 < n < iV. Obviously, it is D„ = and = 1m- These 
discrete symmetries are non-Abelian in the adjoint group space, in the sense that the eigenvalue matrix c"^, 
determined by means of the relation Z)„T"£)„ = c"^T^, is not diagonal. 
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Table 4: Particle spectrum of the inRaton and an SU (2) x -gauged waterfall sectors after inRation. 
The would-be Goldstone bosons of the respective SU(2)x gauge fields are given in the square 
brackets 



Analogously to the U(l)x case, we express the SU(2)x doublets X12 in terms of 
eigenstates of the 2-parities [cf. (4.4)]. In terms of their components, these fields may 
be conveniently expressed as follows: 

±X± = + -L + i±7^^ , (4.24) 

with (X^) — y/2M and {X_) = in the absence of any D-parity violating couphng in the 
theory. Moreover, the SU(2)x -D-terms are given by 

D« = I (xlr^Xi - XJt'^X;^ . (4.25) 

In the D-parity eigenbasis (4.4), they take on the form 

g \ X|r"X_ + X^ r"X+ . for symmetric (a = 1, 3) 
D"=-x<^t ^ (4.26) 

2 [ xIt^X^ + Xlr"X_ , for r" antisymmetric (a = 2) 

Exactly as in the U(l)x case, we find that there are two groups of mass-degenerate fields, 
K- and gi-sector, with masses and rUg given in (4.11) and (4.13), respectively. The 
complete inflaton-waterfall spectrum, along with their Di and D2 parities, is exhibited in 
Table 4. 

The conservation of both 2-parities enforces the stability of all gf-sector particles. 
Instead, if only the L>i-parity, but not £)2, is conserved, then only the i?i-odd particles from 
Table 4 will be stable, and vice versa. Obviously, both Di- and D2-parities need be broken 
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to make all (jf-scctor particles unstable. In Appendix A, we discuss possible mechanisms of 
explicit D-parity breaking for an SU(2)x waterfall-gauge sector. In general, there are two 
mechanisms for breaking D-parity. The first one consists of including higher-order non- 
renormalizable operators in the Kahler potential whose presence explicitly breaks D-parity, 
whilst the second one is very analogous to the U(l)x case. Although a bare FI D-term is 
not possible in non-Abehan theories, effective £)"-tadpole terms may appear after the SSB 
of SU(2)x. The effective iP"-tadpole terms do not break SUSY. They get generated either 
from a non-renormalizable Kahler potential or are induced radiatively, after integrating out 
Planck-scale degrees of freedom. Thus, without excessive tuning, the effective iP"-tadpole 
terms can in general be small of the size required to obtain a second reheat phase in the 
evolution of the Universe. 

Independently of the mechanism which is invoked to break £)-parity, we may in general 
parameterize the g^-sector particle decay rates through the D-parity- violating mass mpi, 
which enters the relation (4.18). In the next section, we will discuss how these relatively 
long-lived g^-sector particles are produced via preheating. 

4.2 Preheating and Thermalization 

After the inflaton field <p passes below a certain critical value 0c ~ M, the so-called waterfall 
mechanism gets triggered. In this case, the inflaton and all other K-sector fields (see 
Tables 3 and 4) oscillate about their true supersymmetric minima: (S) — and (^+) = 
y/2 M. In this waterfall epoch, most of the energy density of the Universe is stored to these 
coherently oscillating ^-sector field condensates and is given initially by — k^M^. During 
the waterfall regime, however, there is an additional mechanism for particle production 
called preheating. 

In general, there are two phenomena associated with the notion of preheating: 

• The first effect of preheating arises from the negative curvature of the potential with 
respect to the K-sector fields. Such a negative curvature corresponds to a negative 
tachyonic mass term in the potential. As a consequence, the particle number within 
infrared modes of momentum less than this tachyonic mass grows exponentially. This 
phenomenon is known as the negative coupling instability or tachyonic preheating [48] . 
Numerical simulations have shown that the field amplitudes suffer strong damping 
during the first oscillation, due to the energy transfer to the infrared modes. In the 
Fo-term hybrid model, only K-sector particles are produced by tachyonic preheating. 
A full study of this process, including thermal equilibration of the K-sector particles, is 
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a highly nontrivial matter and has so far only been achieved for very particular models 
of preheating. Since the fraction of the energy density transferred instantaneously 
to K-sector particles through tachyonic preheating is rather small, compared to their 
initial energy density p^, these model-dependent details fortunately have no dramatic 
impact on the expansion and the thermal history of the Universe. Therefore, we do 
not consider the phenomenon of tachyonic preheating in the Fp-term hybrid model. 

• Particle production may also occur during the coherent oscillation regime, because 
both the K- and ^f-sector particles have masses that can vary very strongly with time. 
This effect is called preheating via a time-varying mass or simply preheating [49,50]. 
As we will show below, a small but significant fraction of the total energy density of 
the Universe p^^ can be transferred, almost instantaneously, to the (yf-sector particles, 
e.g. pg ~ IO^^Pk, for k ~ 10~^. As we illustrated in the beginning of this section 
and will show more explicitly in Section 5, this small fraction of the (^-sector energy 
density is sufficient to alter dramatically the thermal history of the Universe. 

Our interest lies therefore in computing the production energy density pg of the g- 
sector particles via preheating. A key element in such a computation is the profile of the 
time- varying mass of the g^-sector particles, mg{t) = g X+{t) /y/2. The exact time depen- 
dence of mg{t) depends crucially on the dynamics of tachyonic preheating. Comparative 
numerical studies strongly suggest that a sufficiently accurate description of the time evo- 
lution of the g'-sector mass is obtained by [50]^ 



Notice that the time-dependent function mg{t) properly interpolates between the values 
mg{t — > — cxo) = and mg{t — >• oo) = gM that occur in the beginning and the end of the 
waterfall epoch, respectively. 

Given the time-dependent mass (4.27), we may compute the occupation number of 
the fermionic gi-sector modes by solving the Dirac equation 



'''To be specific, the mass-term time- variation stiidied in [50] was for a model with a single field rolling 
from the top of a local maximum of a quartic potential. It was found that a tanh-functional dependence 
accurately captures the evolution of the time- varying mass. Even though the model considered [50] is still 
different from our hybrid inflationary potential, the derived tanh-functional profile for the time-varying 
mass should be regarded as a substantial improvement over the one assumed in [13]. 



mg{t) 



gM r 



tanh(KMi) + 1 



(4.27) 



2 L 



i7°at - 7-k - mg{t) Uh{t) = 0. 
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The solution to the above equation may be expressed by the time-dependent Dirac spinor 
Uh{t) in the chiral representation: 



Uh{t) 




(4.29) 



where is the hehcity two-component eigenspinor for hehcity h — ±. The occupation 
number of Dirac fermions produced via preheating in the true supersymmetric vacuum at 
i — > oo is given by 



nlik) 



1 



1 

+ 2- 



(4.30) 



where k = |k| is the modulus of the 3-momcntum. With the help of (4.30), the /c-mode 
energy density is calculated by p{k) = '^fi'-o^k) n^{k), where u{k) = \Jk'^ + m?g{t oo). 
To obtain a unique solution to the linear differential equation (4.28), we impose initial 
conditions that correspond to a zero occupation number, i.e. n\{k) = 0. These are given 
at t — > —oo by 



uj{k) + hk 



Rh — 



uj{k) — hk 



(4.31) 



2a; V 2a; 

By analogy, the occupation number of the bosonic gf-sector modes are determined by 
solving the Klein-Gordon equation of motion 



0, 



and imposing the initial conditions at t 

1 



-oo. 



'at 



— 1 



(4.32) 



(4.33) 



As in the Dirac case, these initial conditions correspond to vanishing occupation numbers. 
The occupation number of the bosonic modes at f — > oo is given by 

2 



n^(k) 



-u;(k)\^\^ + 



2Lo{k) 



dip 



dt 



1 
2 



(4.34) 



Using the time-dependent mass-term (4.27), along with the initial conditions (4.31) 
and (4.34), one may obtain analytical expressions in terms of hypergeometric functions [50], 
for the particle production between t — > — oo and t — > oo. For k <^ g and k <^ gM, these 
analytical expressions reduce to 

2 



n{k) 



(4.35) 
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where the sign + apphes for n{k) = n^ik) and the sign — for n{k) = n^{k). Recalhng that 
there are 2 hehcity states for a (yf-sector fermion and 4 real degrees of freedom for a ^f-sector 
boson, we may calculate the occupation number of all ^f-sector modes as follows: 

ng{k) = iVb (^Yl + j ' (4-36) 

where A^t, is the number of broken generators of the waterfall gauge symmetry. In particular, 
it is A^b = 1 for U(l)x and A^b = 3 for SU(2)x [cf. Tables 3 and 4]. Since the produced 
particles are non-relativistic, i.e. k <^ gM, their occupation number distribution ng{k) can 
easily be integrated to give the total energy density carried by the (^-sector fields, i.e. 



Pg ^ gM 



oo 

k^dkugik) ^ 2.1 X 10-2 iVb . (4.37) 



Pk Pk 27r2 



Here = k^M^ is the energy density of the K-sector particles shortly before the waterfall 
transition. Equation (4.37) will be a valuable input for the next section to compute the 
true reheat temperature Treh of the Universe, which arises from the combined effect of the 
K- and gf-sector particle decays. 



5 Coupled Reheating and Gravitino Abundance 

In the previous section, we have seen that the g^-sector particles, e.g. ipg, R- and V^, can 
be abundantly produced during the preheating epoch. Assuming that they dominate the 
Universe at some later time, their decays induced by the small D-parity violating couplings 
will give rise to a second reheat temperature, which we denote here by Tg. As we will 
show in this section, the large entropy, which is released by the late decays of the gi-sector 
particles, will be sufficient to dilute the gravitinos to levels compatible with BBN fimits. 

More expUcitly, we present a detailed numerical analysis of the gravitino abun- 
dance Yg, where the combined effect of the k- and gi-sector particle decays is carefully 
taken into account. As we mentioned in the Introduction, we call such a two-states' mech- 
anism of reheating the Universe coupled reheating. In Section 5.1, we set the BEs relevant 
to coupled reheating and give numerical estimates of the gravitino abundance Yq and the 
energy densities p^, pg and prad related to the k- and g^-sector particles and their radiation, 
respectively. In Section 5.2, we present a semi-analytic approach to BEs, where useful 
approximate expressions for Yq are obtained. Finally, in Section 5.3 we derive gravitino 
abundance constraints on the theoretical parameters. 
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5.1 Boltzmann Equations 



The number density Uq of gravitinos, the energy density (pg) of the k (gi)-sector particles 
and the energy density prad of the radiation produced by their decays satisfy the following 
system of BEs [53]: 

fiQ + SHuQ = CqT^, 
Pk + 3Hp^ = -TkPk, 

Pg + SHpg = -^gPg, 
Prad + 4-f/"Prad = ^nPn + ^ g Pg , (5-1) 

where a dot on ng, p^.g and prad indicates differentiation with respect to the cosmic time t. 
The quantity Cq{T) is a collision term for gravitino production calculated in [51,52] and 
the Hubble parameter H is given by 



1 / xl/2 

H = —= mgng + p« + Pg + Prad , (5.2) 

A/.i77?,oi V / 



where is the mass of the gravitino G. In addition, the temperature T and the entropy 
density s may be determined through the relations: 

Prad = —g.T\ s ^ —g.T\ (5.3) 

where g^{T) is the effective number of degrees of freedom at temperature T. Since the 
initial temperature is Tin ^ z^^, it is = 240 for all T > Msusy- 

Here we should note that in BEs (5.1) we have neglected the collision terms related 
to the self-annihilation of (/-sector particles. Their thermally averaged cross section times 
velocity, {aannv), is estimated to be 

(cXann^^) < 10"^^ GeV-% (5.4) 

which is numerically negligible. 

The numerical analysis of the BEs (5.1) gets simphfied by absorbing the Hubble 
expansion terms into new variables. To this end, we define the following dimensionless 
quantities [54]: 

fd = "-G^^ /« = fg = /rad = PradO^- (5.5) 

where a is the usual expansion scale factor of the Universe. We also convert the time 
derivatives to derivatives with respect to the logarithmic time ln(a/ai) [55], where oi is 
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some initial or reference value for the scale factor a. With the above substitutions, the 
BEs (5.1) may be re-written as 

Hf~ = Cg T^a^ , 
Hf'n = - Tk/k , 
Hfg = -^gfg, 
HfL = r^/<^« + r,/,a, (5.6) 

where the prime now denotes differentiation with respect to In (a/aj). Correspondingly, 
the Hubble parameter H and temperature T may now be expressed in terms of the newly 
introduced variables (5.5) as follows: 

" = ('"«^« + ^ ^ ""^"^y"' ^ = i'^T ■ 

The transformed system of BEs (5.6) can be numerically solved by imposing the following 
initial conditions: 

Ujaf = K^M\ fgjaf = 2.1x10-^ gK'M\ /rad,i = , (5.8) 

where the subscript I refers to quantities defined at ln(a/ai) = 0. Notice that the initial 
value fgj af is equal to the energy density p^j of the gf-sector particles produced during 
preheating and is given in (4.37) . 

In Fig. 7(a), we present numerical estimates of the cosmological evolution of energy 
densities PK,g,ra.A as functions of the temperature T in a double x-y logarithmic plot, where 
is represented by a dark grey line, pg by a grey line and piad by a light grey line. As an 
example, we use M = 0.7 x 10^^ GeV, p = A = k = 10"^ and mpi/M = 4.3 x 10"'^ (bold 
lines) and m-pi/M = 10~^ (thin lines). Since prad is affected very little for the larger value 
of mpi/M, it has not been added to the plot. The intersection point of the T-dcpcndent 
functions p^ and p^ad signals the completion of the /t-sector particle decays. For the specific 
example, this point occurs at = 3.2 x lO"*^^ GeV. For nipi/M = 4.3 x 10^^ GeV, we obtain 
two more intersections: one for T = Tgq ~ 3.9 x 10^ GeV where pg{Tcq) = Prad(^eq) and 
another one for T = Tg c:^ 200 GeV, where the (y'-sector particles have practically decayed 
away and Pg{Tg) = Pra.d{Tg). Thanks to the huge entropy release in this case, the gravitino 
abundance = n^/s gets sharply decreased from about 2.2 x lO^"*^^ to 2.4 x lO^"*^^. This 
dramatic reduction of is shown in Fig. 7(b). On the contrary, if m-p\/M = 10~^, no 
intersection of Pg with p^ad takes place and, in consequence, no phase of second reheating 
occurs. This is also illustrated in Fig. 7(a), where the dependence of pg is displayed by a 
thin line. As can be seen from Fig. 7(b), the gravitino abundance remains unsuppressed 
in this case, i.e. Yq ~ 10"^'^. As we will see below in Section 5.3, such large values of Yq 
are in gross conflict with BBN constraints. 
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Figure 7: The evolution as a function of log T of the quantities: (a) log pi with i = k (dark 
grey Une), i = g (grey Une), i = rad (Ught grey Une) (b) G yield, Yq. In both cases, we take 
M = 0.7 X 10^^ GcV, p = \ = K = 0.001 and mpi/M = 4.3 x 10"^ GeV (bold lines) and 
mpi/M = 1 X 10-2 (thin lines). 
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5.2 Semi-analytic Approach 



We now present a more intuitive and rather accurate approach to the dynamics of coupled 
reheating, and find approximate analytical expressions that describe the evolution of the 
energy densities p K,g, raA- In addition, we derive the conditions that ensure the existence of 
a second reheat phase in the evolution of the Universe. Finally, we estimate the gravitino 
abundance Yq due to coupled reheating. 

Shortly after infiation ends, the energy of our observable Universe is dominated by 
the infiaton S and the other re-sector particles, with an initial energy density p^j = k^M^. 
As we schematically illustrated in Section 4.2, the re-sector particles decay into relativistic 
degrees of freedom, producing an energy density prad- The energy density Pg of the gi-sector 
particles is subdominant at these early stages after the first reheating due to the re-sector 
particle decays. In fact, for temperatures T > Teq, where Teq is the first intersection point 
of the T-dependent functions prad and pg [see (5.10)], the evolution of all relevant energy 
densities may be approximately described as follows: 



Pk = Pk,i (a/ai) , pg = pg^i (a/ai) , p^ad = Prad{T^) (T/T^) . (5.9) 

As mentioned above, the T-dependent function prad may first cross the corresponding pg 
at T — Teq, where 

Prad(7;q) = P,(req) ■ (5.10) 

Using the fact that PradiT^) = Pk(7k) and assuming that the Universe expands isentropically 
with a oc T"^ when Teq <T< T^, we obtain from (5.10) the approximate relation 

Teq ^ r« ^ . (5.11) 

In deriving (5.11), we have also made use of (5.9). 

A second reheat phase in the evolution of the Universe takes place, only if Tg < Teq, 
where Tg is the naive reheat temperature due to the (^-sector particles decays [see (5.20) 
below]. To better explore the consequences of this last condition, wc use the abbreviation 
g-DAD ((yf-DBD) to indicate whether the (^-sector particles Decay After (Before) the Dom- 
ination of their energy density. With the aid of (5.11), the following two conditions for 
g-DAD and gf-DBD may be deduced: 

^ < ^ (5-DAD), ^ > ^ (5-DBD) . (5.12) 

These two possible scenarios are illustrated in Fig. 7 for mpi/M — 4.3 x 10"*^ (mpi/M — 
10~^), where the bold (thin) lines correspond to g-DAD (g-DBD). 



43 



The gravitino abundance Yq can be calculated by simply integrating /- that occurs 
in the first BE of (5.6) and using the fact that = f^/sa^. It turns out that the main 
contribution to comes from the integration after the commencement of the radiation 
dominated era, i.e. for T < T^. The so-derived formula reproduces rather accurately the 
one presented in [51] in the massless gluino limit, where 

Note that (5.13) is only valid for the (yf-DBD case. 

The situation is different for the g-DKD case, where a drastic reduction of the gravitino 
abundance, caused by the huge entropy release from the (yf-sector particle decays, takes 
place. In this case, the gravitino abundance Y~ may be estimated in the following way. We 
first notice that ^ 

Y9 ^ YH ^(-^eq) Q (^eq) z^- 

o G s{T,)a^{T,) ■ ^ 
Then, with the help of (5.3) and (5.9), we may obtain the relation 

g(Teq)a3(Teq) ^ f T^V f ^ ^ 

s{T,)a^{T,) \Tj \p,{T,)) ' ^ " ^ 

Substituting the respective expressions of (5.15), (5.13) and (5.11) into (5.14) yields 

where we have used (4.37) to derive the last approximate equality. We have checked that 
the semi-analytic formula (5.16) is in remarkable agreement with numerical estimates in 
the g-DAD regime. 

Finally, we should comment on the fact that the number of e-folds, gets modified 
in the g-DAD case, because of the occurrence of a (?-sector-mattcr dominated era. Making 
use of standard methods [6,8], we are able to determine A4 at the WMAP pivotal point 
ko = 0.002 Mpc^^ by the following relation: 

K = 22.6 + hni^M'') + llaTg + ^ In ^ . (5.17) 
6 3 3 pg^i 

This result, however, does not crucially alter the value of A^e, which remains close to 55 — 60 
in the g-DAD well. Interestingly enough, Y~ and J\fe do not directly depend on 

given in (4.15). In fact, in the g-DAD case, Y~ and are fully independent of the 
supcrpotential couplings A and p, and only have a mild linear and logarithmic dependence 
on K, respectively. As we will discuss below, it is this last property that leads to a significant 
relaxation of the strict gravitino constraints on these couplings, when compared to the g- 
DBD case. 
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5.3 Gravitino Abundance Constraints 



In order to avoid destroying the apparent success between the standard theory for BBN 
and observation, gravitinos must have an abundance Yq below certain upper hmits, which 
crucially depend on their decay properties [51,56]. Some representative upper bounds 
on Yg, obtained in a very recent analysis [51], are 

f 10-15^ for mg~ 360 GeV, 
10-l^ for mp,~ 600 GeV, 

5.18 

10-^^ for ruQ ~ 7.5 TeV, ^ ' 

10-^2^ for m^~9.3TeV. 
The above bounds pertain to the less restrictive case of a gravitino that decays with a small 
branching ratio Bh — 0.001 into hadronic modes. For the g^-DBD case discussed above, the 
upper limits (5.18) imply the corresponding stringent bounds on Treh: 

9 X 10^ GeV, for tuq ~ 360 GeV , 
7x10^ GeV, for rriQ ~ 600 GeV , 
T.eh ^ < 7 X GeV, for tuq ~ 7.5 TeV , ^^'^^^ 
7 X lO'^ GeV, for ~ 9.3 TeV . 

The aforementioned upper limits lead to serious constraints on the basic couplings 
A and p, usually forcing them to acquire very small values, i.e. k. X, p ^ 10~^. For the 
standard F-term hybrid model within mSUGRA and with a soft SUSY-breaking tadpole 
parameter a^ = 1 TeV, the requirement of accounting for the observed power spectrum P-n, 
with a number of e-folds jVe = 50-60, implies that k > 10'^ and T« = T^eh ^ 9 x 10^ GeV. 
Such a high lower bound on Treh invalidates all the limits presented in (5.19), thereby ruling 
out the above F-term hybrid model. 

The above situation, however, changes drastically in the Fp-term hybrid model with 
small £)-parity violation, e.g. due to the presence of a subdominant FI D-term. This 
corresponds to the g^-DAD case described in the previous subsection, where the upper 
bounds (5.18) translate, by means of (5.16), into upper bounds on mpi/M for k > 8 x 10~^. 
The required size of the £)-parity violating parameter mpi may naively be estimated using 
a relation very analogous to (4.15), viz. 

90 ^ 



T, = [^J , (5.20) 

where is the decay width of a g'-sector particle and is given in (4.18). If we solve (5.20) 
for the ratio mpi/M, we obtain 

« 8.4 . 10- X f " (^) (15^) (5.21) 

M \ g J \lO^ GeV J \ M J ^ ^ 



45 



T (GeV) 



10' 10" 10' 10' 10' 10' 10' 10'° 10'' 10'' 10'' 




ig f ' ' ' ' ' 

10"' 10"' 10"* 10"' 10"' 



rripi / M 

Figure 8: The dependence of log on mpi/M, for k = 10~^ (dark grey line), k = 10~^ (grey- 
line) and K = 10~^ (light grey line). 



For second reheat temperatures Tg of cosmological interest, i.e. 0.2 TeV ^ Tg < 10^ GeV, 
the following double inequahty for M — 10^^ GeV may be derived: 

4 X 10-^ < ^ < 10-3 . (5.22) 

The lower bounds on Tg and m-pi/M result from the requirement that thermal electroweak- 
scale resonant leptogenesis be successfully realized. More details are given in Section 6. 

A numerical analysis of the gravitino abundance predictions and constraints related to 
the ^f-DAD scenario has been performed in Figs. 8 and 9, respectively. Our numerical results 
apply equally well to both mSUGRA and nmSUGRA scenarios. In detail. Fig. 8 shows 
logy^ as a function of mpi/M and Tg, for the different values oi k, — 10~^, 10"^ , 10"^, 
while M is fixed by the usual inflationary constraints on and P-jz, ior k — \ — p 
and ch — ^- The different hues in Fig. 8 terminate at high values of mpi/M, since the 
inequality Tg < T^ does no longer hold. The lowest value of mpi/M is determined by the 
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10"* 10"' 10"' 

K 

Figure 9: The allowed region on the (mpi/M, k) plane for Yg < 10~^^ (black area), < 10~^^ 
(light grey area), < 10~^^ (grey area) and Y^ < 10'^^ (dark grey area). 

condition Tg > 200 GeV, which results from the aforementioned requirement that thermal 
electroweak-scale resonant leptogenesis is successfully realized [16,19,20]. 

In Fig. 8, we also observe the two regimes: g-DBD and g-DAD. In the g-DBD regime, 
logFg remains constant for given k up to some value m-pi/M. For example, for k = 10~^, 
logFg is constant for mpi/M > 10^"^. This result is consistent with (5.13). For smaller 
values of m-pi/M, one enters the g-DAT> regime. In this case, logFg decreases rapidly, 
as rriYi/M, or equivalent ly Tg, decreases. This behaviour of Yq is expected on account 
of (5.16). Also, in agreement with (5.16), the reduction of Yq becomes more drastic for 
larger values of k. 

In Fig. 9 we delineate the allowed regions on the {K,m-pi/ M) plane for the discrete 
values of Yq = 10"^^, 10"^"^, 10~^^, 10~^^, for k > 8 x 10~^. The upper boundaries of the 
various areas are obtained using (5.16). For k < 8x 10~^, we are in the g-DBD region, where 
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we obtain 10~^^ < < 10""^^, almost independently of m-pi/M [cf. (5.13)]. Therefore, we 
only display values for itlfi/M, for which g-DAD becomes relevant. We observe that the 
most stringent limit on can still be fulfilled for k > 10~^ and m-pi/M < 10^®. Such 
large values of k would have been excluded from naive estimates of the K-sector reheat 
temperature due to the K-sector particle decays. According to our analysis in this 
section, however, these large values of k, A and p of order 10~^-10~^ are allowed within the 
Fo-term hybrid inflationary model. As we will see in the next section, this is a distinctive 
feature of the Fp-term hybrid model that opens up novel possibilities in solving the CDM 
problem. 

At the end of this section, we wish to comment on a possible F^i-term hybrid scenario, 
where the K-sector particles can decay directly into the g^-sector ones. This can happen, for 
example, if m«; > 2mg or equivalently when n > \/2g. Since the gauge coupling g of the 
waterfall sector must be smaller than 0.1 in this case, it would be difficult to embed such 
a Fp-term hybrid scenario into a GUT. The energy density transferred from the K-sector 
particles into the g^-sector ones may be calculated by 

^ = -|- B.^, . (5.23) 

Here B^^g denotes the branching ratio of the decays of the k- to (^-sector particles. Assum- 
ing conservatively that -B^^g ~ 10^^ and k = 2g, we obtain an estimate for the gravitino 
abundance ~ 10~^® for m-pi/M < 10~^, thereby rendering gravitinos quite harmless. 



6 Baryon Asymmetry and Cold Dark Matter 

In this section we briefly discuss further cosmological implications of the Fp-term hybrid 
model for the BAU and the CDM. 



6.1 Resonant Flavour-Leptogenesis at the Electroweak Scale 

Earlier studies of the BAU in supersymmetric models of hybrid inflation have mainly been 
focused on scenarios of non-thermal leptogenesis [57] , with an hierarchical heavy Majorana 
neutrino spectrum, e.g. ruNi < mjVa <^ niNa- The simplest model of this type is obtained 
by identifying the waterfall gauge group with \J{1)b-l, which allows the presence of the 
operator X2X2NiNj/m-p\ in the superpotential. Notice that such a term is forbidden 
in the F^-term hybrid model by virtue of the it! symmetry (2.6). In the non-thermal 
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leptogenesis model, the reheat temperature consistent with the observed BAU tje = 6.1 x 
10~^° and low-energy neutrino data is estimated to be [12] 

where the superpotential couplings A, p are set to zero, li X — k, the lower bound (6.1) gets 
larger roughly by a factor 20. It is obvious that in this generic non-thermal leptogenesis 
scenario, the gravitino constraint on Treh favours rather small values of k and A, e.g. k, X ^ 
10~^ for Treh ~ 10^ GeV. As was discussed in Section 3.2, however, such small values of k, 
introduce strong tuning at a less than 1% level to the horizon exit values of the inflaton 
field 0exit in a nmSUGRA scenario that accounts for the recently observed value of the 
spectral index given by (3.10). Moreover, the success of this scenario rehes heavily on 
the assumption that there is no other source of baryogenesis, e.g. through the Affleck- 
Dine mechanism, nor of entropy release, e.g. from possible late decays of moduli or flaton 
fields [58], between the energy scales mjvi (^ Treh) and the electroweak phase transition. 

In the Fp-term hybrid model, non-thermal leptogenesis is not possible for one of the 
reasons mentioned above. The late decays of the g'-sector (D-odd) particles generally lead 
to an enormous entropy release, so that not only gravitinos, but also any initial lepton- 
number excess will be diluted to unobservable values. However, as has already been dis- 
cussed in [13], the F^-term model can realize electroweak-scale resonant leptogenesis [16], 
if the couphng of the inflaton superfield S to the respective right-handed neutrinos Ni is 
very close to an S0(3)-symmetric form, i.e. pij ~ p I3. This will give rise to 3 nearly heavy 
Major ana neutrinos of mass mjv and so would enable a successful realization of the reso- 
nant leptogenesis mechanism at the electroweak scale. The required S0(3)-breaking may, 
for example, originate from renormalization-group (RG) [59] or possible GUT threshold 
effects [16,60]. 

An order of magnitude estimate of the final BAU rjB, including single lepton flavour 
effects, may be obtained as [16, 20] 

rjB ^ - 10-' X r{Tjm^) J] J] S'^^ , (6.2) 

where 

r,i _ r(Ar, ^ Lz$) + T{N, ^ Lf $t) 

" H{T = m^) ^^■'^^ 

is a lepton-flavour dependent wash-out factor, which quantifies in a way the degree of in- 
or out-of-equilibrium of the decay rates of the heavy Majorana neutrino mass eigenstates 
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Figure 10: Numerical estimates of the BAU for a scenario with = 250 GeV and for different 
initial lepton- and baryon-number abundances, rj^ and rf^, assuming an initial thermal distribu- 
tion for the heavy Majorana neutrinos, i.e. r/]^^ 23 ~ ^' "^^^ horizontal grey line shows the BAU 
needed to agree with today's observed value. 

Ni {i = 1,2,3) into the SM-like Higgs doublet $ and the lepton doublet L/ (/ = e, /x,r). 
The remaining i^T-factors in (6.2) are defined with the help of K'^^, as follows: 

3 

1=1 Ni 

The parameters 5jy. denote the different lepton-flavour asymmetries related to the decays 
Ni ^ Li^ and are defined by 



Finally, the prefactor r{Tg/mN) in (6.2) takes care of a possible dilution effect on the BAU 
that might be caused by the entropy release of late ^f-sector particle decays. This dilution 
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effect is only relevant, if the second reheat temperature Tg is smaller than the leptogenesis 
scale rriN- Employing standard arguments of thermodynamics, one may estimate that 



Instead, if Tg 3> m^v, the dilution factor r{Tg/m.N) approaches 1 and can therefore be 
omitted. 

In Fig. 10, we display numerical estimates of the BAU for a resonant leptogenesis 
scenario with ttin — 250 GeV and an inverted hierarchical light- neutrino spectrum. For 
a detailed discussion of the heavy and light neutrino spectra of this model, the reader 
is referred to [16]. As can be seen from Fig. 10, one advantageous feature of resonant 
leptogenesis is that the final baryon asymmetry tib does not sensitively depend on any pre- 
existing lepton- or baryon-number abundance, 77™^ or 77^. For instance, assuming an initial 
thermal distribution for the heavy Majorana neutrinos, i.e. riNi2 3 ~ primordial 
baryon asymmetries rj^ ^ 10'^, we observe that the final 77B is practically independent 
of the initial conditions, once the relevant particle-physics model parameters, such as the 
heavy Majorana masses and their respective Yukawa couplings, are fixed. 

It is important to comment here on the fact that the above property of the indepen- 
dence of the BAU on the initial conditions does not necessarily get spoiled, if the second re- 
heat temperature Tg happens to be smaller than the resonant leptogenesis scale uin- In this 
case, one only needs to make sure that the entropy dilution suppression factor ~ (Tg/mjv)^ 
does not lead to a significant reduction of the BAU. Therefore, we have rather conserva- 
tively assumed throughout our numerical analysis in Section 5 that ^ mjv ~ 250 GeV, 
even though Tg could still be somewhat smaller than the resonant leptogenesis scale mjv. 

Another point that deserves to be clarified here is the physical significance of lepton- 
fiavour effects on the BAU. In general, there are two sources of lepton fiavour: (i) the 
charged lepton Yukawa couplings hi and (ii) the neutrino Yukawa couplings h^j . The former 
has been extensively discussed in the literature [61] and may affect the predictions for the 
BAU by up to one order of magnitude, depending on the scale of leptogenesis. For our 
electroweak-scale leptogenesis scenario, these effects are not significant, since all charged 
lepton Yukawa couplings mediate interactions that are in thermal equilibrium. The second 
source of fiavour effects is due to neutrino- Yukawa couplings h^j and has been studied 
only very recently in [16,20,62]. The effect on the BAU is most relevant when the heavy 
Majorana neutrinos get closer in mass. In models of resonant leptogenesis, neutrino- Yukawa 
coupling effects can have a dramatic impact on the predictions for the BAU, enhancing its 
value by many orders of magnitude [16,20]. 




(6.6) 
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This last fact opens up new vistas in the model-building of scenarios that can be 
phenomenologically more accessible to laboratory experiments. For instance, if a certain 
hierarchy among the Yukawa-neutrino couplings h\j is assumed, e.g. h^2 = ^K'i ~ ^0^^ ^ 
and /ij'j = 10~®-10~^ ~ he, resulting from the approximate breaking of some global 11(1); 
symmetry, the required BAU can still be generated successfully from an individual lepton 
number asymmetry, namely L-,- in this case. For this particular model of resonant r- 
leptogenesis, the values of the X-factors defined in (6.3) are: 

x;,^^3 ~ 10, ~ 30, x^j;^ ~ 10^°. (6.7) 

Given that the leptonic asymmetry is 5]^^ ~ 10~^, one can estimate from (6.2) that the 
right amount of baryon asymmetry is produced, with r]B ~ 10~^. This is also shown in 
Fig. 10. Instead, older approaches to BEs that do not appropriately treat lepton flavour 
effects via the neutrino Yukawa couplings h'^j would have predicted a value that would have 
been short of a huge factor ~ 10~^ [16,20]. 

As can be seen from the above example, the lepton-flavour directions Lg,^ orthogonal 
to Lr can involve large neutrino Yukawa couplings of order 10^^. Such couplings can 
give rise to distinctive signatures in the production and decay of electroweak-scale heavy 
Majorana neutrinos at high-energy colliders, such as the LHC [63] , the International Linear 
e+e" Collider (ILC) [64] and other future colliders [65]. Moreover, electroweak-scale heavy 
Majorana neutrinos can give rise to phenomena of lepton flavour and/or number violation, 
such as the neutrinoless double-beta decay (Oz//3/5), the decays — > 67 [66], /i —>■ eee, /i — > e 
conversion in nuclei etc [67-70]. A detailed discussion of the low-energy phenomenology of 
resonant leptogenesis models may be found in [16]. 

6.2 Thermal Right- Handed Sneutrinos as CDM 

An interesting feature of the Fp-term hybrid model is that it!-parity is conserved, even 
though the lepton number L, as well as B — L, are explicitly broken by the Majorana 
operator | pSNiNi. In fact, in our model, all superpotential couplings either conserve the 
B — L number or break it by even number of units. For example, the coupling p breaks 
explicitly L, along with B — L, by 2 units. Since the i?-parity of each superpotential 
operator is determined to be = (^—I^H^-l) _ _|_-|^^ ^j^g Fp-term hybrid model conserves 
i?-parity. As a consequence, the LSP of the spectrum is stable and so becomes a viable 
candidate to address the CDM problem of the Universe. 

In addition to the standard CDM candidates of the MSSM, e.g. a stable neutralino, 
it would be interesting to explore whether thermal right-handed sneutrinos as LSPs could 
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solve the CDM problem. Before we estimate their relic abundance, we first observe that 
light right-handed sneutrinos may easily appear in the spectrum. Ignoring the small 
neutrino- Yukawa coupling terms, the right-handed sneutrino mass matrix is written 
down in the weak basis (iV"i^2,3, -^1,2,3) • 

2 _ f / p'vl + Ml pApVs + pXvuVd \ , . 

2[ pA;vs + pXv^va p'vl + Ml ) ' ' 

where Vs = (S), Vu,d = {Hu,d) and M"- is the soft SUSY-breaking mass parameters associ- 
ated with the sneutrino fields. The sneutrino spectrum will then consist of 3 heavy (light) 
right-handed sneutrinos of mass 

p'^vl + M| + (-) (^pApVs + pXvuVd^ . 

Hence, the 3 light sneutrinos can act as LSPs, which we collectively denote them by Nlsp- 

Recently, the possibility that right-handed sneutrinos are the CDM was considered 
in [71]. This recent analysis showed that thermal right-handed sneutrinos have rather 
high relic abundances and will generally overdose the Universe in a supersymmetric exten- 
sion of the MSSM with right-handed neutrino superfields Ni and bare Majorana masses 
{mM)ijNiNj. The underlying reason is that because of the small Yukawa- neutrino couplings 
/if, , the self- and co- annihilation interactions of the sneutrino LSP with itself and other 
MSSM particles are rather weak. These weak processes do not allow the sneutrino LSP 
to stay long enough in thermal equilibrium before its freeze-out temperature, such that its 
number density gets reduced to a level compatible with the CMB data, i.e. Qdm^^ ~ 0.15. 
Instead, the predicted values turn out to be many orders of magnitude larger than 1. 

In the Fo-term hybrid model, however, there is a new interaction that can make the 
right-handed sneutrinos annihilate more efficiently. This is the quartic coupling ^ 

A'if = \xpN;N;H^Hd + h.c. (6.9) 

It results from the F-term of the inflaton field: Fs ~ \ pNiNi + XHuHd- To assess 
the significance of the interaction (6.9), we estimate the relic density of A^lsp in different 
kinematic regions. 

^The implications of a generic quartic coupling of the same form for the CDM abundance and detection 
was studied earlier in [72,73] within the context of a simple non-SUSY model. These studies will not be 
directly applicable to our more elaborate case of a supersymmetric scenario with right-handed sneutrinos. 
However, we have used their results to check our qualitative estimates for the CDM abundance. 
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We first consider the self-annihilation off-resonant process A^lsp^lsp ~^ {Hy) — >■ 
VT+iy", which occurs when ?^jvlsp ^ ^w- A simple estimate yields 

ikuh ~ (^^j (^M^j ■ (6.10) 

To obtain an acceptable relic density, we need relatively large p and A couplings, i.e. p, A ^ 
0.1^. Such values go in opposite direction with those obtained by requiring successful 
inflation with a red-tilted spectrum. Therefore, as far as inflation is concerned, they signify 
the necessity of going well beyond the minimal Kahler potential. 

The above situation may slightly improve for ^tvlsp ^ ^w, in large tan/3 scenarios 
with light Higgs bosons that couple appreciably to 6-quarks [74]. In particular, in the 
kinematic region Mh^ ~ 2m^^^p, the self-annihilation process Ni^spNhSP {Hu) — >■ hb 
becomes resonant, and the above estimate modifies to 

Qdm/i' ~ 10-4xS^i(if,-.iVLSpiVLSp) X (y^^^) ■ (6.11) 

Consequently, if the couphngs A, p are not too small, e.g. A, p > 10~^, the LSP right-handed 
sneutrinos A^lsp can efficiently annihilate via a Higgs resonance into pairs of &-quarks 
in this kinematic region, thus obtaining a relic density compatible with the CMB data. 
A detailed study of the thermal right-handed sneutrino as CDM could be given elsewhere. 



7 Conclusions 

We have analyzed the cosmological implications of a novel F-term hybrid infiationary 
model, in which the inflaton and the gauged waterfall sectors respect an approximate 
discrete symmetry which we called here D-parity. The approximate breaking of D-parity 
occurs explicitly either through the presence of a subdominant FI D-term or through non- 
renormalizable operators in the Kahler potential. For brevity, this scenario of infiation was 
termed F/p-term hybrid infiation. One of the most interesting features of the model is that 
the VEV of the infiaton field closely relates the p-parameter of the MSSM to an SO (3) 
symmetric Majorana mass ttin- If A ~ p, this implies that p ~ mAr, so the Fci-term hybrid 
model may naturally predict lepton-number violation at the electrowcak scale. 

^An upper bound on the product pA, although somewhat model-dependent, can be derived from exper- 
imental limits on the flux of energetic upward muons that occur in the possible detection of CDM using 
neutrino telescopes [73]. Our initial estimates indicate that it should be pX < 0.03 for m^^^^ ~ 50 GeV, 
which is not a very rectrictive bound. 
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Before summarizing the cosmological and particle-physics imphcations of the Fp-term 
hybrid model, it might be interesting to list our basic assumptions pertinent to inflation 
and to the model itself: 

(i) The standard assumption for successful hybrid inflation is that the inflaton fleld 
should be displaced from its true minimum at the start of inflation, whereas all other 
scalar flelds in the spectrum must have zero VEVs [c.f. (2.14)]. In a nmSUGRA 
scenario of hybrid inflation, however, additional tuning is required beyond the above 
standard assumption. The horizon exit values of the inflaton fleld 0exit have to be close 
to the value 0max, at which the inflationary potential has a maximum. Nevertheless, 
such a tuning is not so strong, i.e. (^max — 0exit)/0exit ~ 10%, as long as k ^ 10~^. 

(ii) Although there may exist several ways of breaking D-parity explicitly, we have con- 
sidered here two possibilities to motivate the required small amount of /^-parity 
violation. As discussed in Appendix A, we first considered the case where D-parity is 
broken by a subdominant FI Z^-term, which is induced radiatively after heavy degrees 
of freedom have been integrated out. Another minimal way would be to introduce 
non-renormalizable operators in the Kahler potential that break £)-parity explicitly. 

(iii) In order to be able to realize thermal resonant leptogenesis at a low scale, the coupling 
matrix p^j is assumed to be close to S0(3) symmetric, i.e. pij ~ pl^. 

The Fo-term hybrid model has several cosmological implications that may be sum- 
marized as follows: 

• The model can accommodate the currently favoured strong red-tilted spectrum with 
ng — 1 ~ —0.05 [4,5], if the radiative corrections dominate the slope of the inflationary 
potential and a ncxt-to-minimal Kahler potential is assumed, where the parameter ch 
is in the range 0.05-0.2. The radiative corrections dominate the slope of the potential, 
if the superpotential couplings, k,, A, p, lie in a certain interval: 10~^ < k. A, p < 10~^. 
In addition, the actual value of the power spectrum Pn and the required number of 
e-folds, N'e ~ 55, provide further constraints on these couplings and the SSB scale 
M of the waterfall gauge symmetry. For example, for M ^ 10^^ GeV, one finds the 
allowed parameter space: k < A, p < An, for k, ~ 10"'^-10~^ and 0.05 < c^f < 0.1. 

• For Fo-term hybrid models with spontaneously broken U(l)x gauge symmetry, the 
non-observation of a cosmic string contribution to the power spectrum at the 10% 
level implies an upper bound on the superpotential coupling k, i.e. k < 10~^. This 
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strict upper bound on k can be weakened by one order of magnitude in a nmSUGRA 
model of Fp-term hybrid inflation, with k, = \ = p and ch = 0.14. On the other hand, 
this upper limit can be completely evaded, if the watefall sector of the F^i-term hybrid 
model realizes an SU(2)x local symmetry that breaks completely to the identity I, 
i.e. SU(2)x — > I. In this case, not only cosmic strings but any other topological 
defects can be avoided, such as monopoles and textures. As we outlined in Section 2, 
GUTs, such as those based on the exceptional groups E(6) and E(7), have breaking 
patterns that contain SU(2)x subgroups uncharged under the SM gauge group and so 
are able to realize F^-term hybrid inflation devoid of monopoles and cosmic strings. 

• To avoid overproduction of gravitinos, one needs to impose a strict upper limit on the 
reheat temperature Tj-ch obtained from the perturbative infiaton decays, i.e. Tj-ch ~ 
10^°-10^ GeV. This upper bound depends on the decay properties of the gravitino 
and gives rise to tight constraints on the basic theoretical parameters k,, A and p, 
i.e. K,, X, p < 10^^. However, these tight limits may be significantly relaxed by 
considering the late decays of the so-called (^-sector particles which are induced by 
small /^-parity violating couplings that may result from either a subdominant FI 
D-term or non-renormalizable Kahler potential terms. These (yf-sector particles are 
produced during the preheating epoch, and if they are abundant, they will lead to a 
second reheating phase in the evolution of the early Universe, giving rise to a rather 
low reheat temperature, even as low as 0.3 TeV. In this case, the enormous entropy 
release from the g^-sector particles may reduce the gravitino abundance below the 
BBN hmits discussed in Section 5. 

• After the infiaton S receives a VEV, one ends up with 3 nearly degenerate heavy 
Majorana neutrinos with masses at the electroweak scale. As we discussed in Sec- 
tion 6, this opens up the possibility to successfully address the BAU within the 
thermal electroweak-scale resonant leptogenesis framework, in a way independent of 
any pre-existing lepton- or baryon-number abundance. 

• The Fo-term hybrid model conserves i?-parity, in spite of the fact that the lepton 
number is explicitly broken by the Majorana operator | p SN^Ni. This is so, because 
all superpotential couplings either conserve the B — L number or break it by even 
number of units. The aforementioned Majorana operator breaks explicitly L, as well 
as S — L, by 2 units. Consequently, the LSP of the spectrum is stable and so qualifies 
as candidate to address the CDM problem. The new aspect of the Fo-term hybrid 
model is that thermal right-handed sneutrinos emerge as new candidates to solve this 
problem, by virtue of the quartic couphng: | XpN*N*HuHd + H.c. This new quartic 
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coupling results in the Higgs potential from the F-terms of the inflaton field, and it is 
not present in the more often-discussed extension of the MSSM, where right-handed 
neutrino superfields have bare Majorana masses. Provided that the couplings A and 
p are not too small, e.g. A, p > 10~^, the LSP right-handed sneutrinos -/Vlsp can 
efficiently annihilate via a Higgs resonance Hd into pairs of 6-quarks, in the kinematic 
region Mh^ ~ 2m ^^^^ , and so drastically reduce its relic density to values compatible 
with the CMB data. 

In addition to the above cosmological implications, the Fp-term hybrid model has 
a rich particle-physics phenomenology. The main phenomenological characteristics of the 
model are: 

(a) It is straightforward to embed the F^-term hybrid model into minimal or next-to- 
minimal SUGRA, where the soft SUSY- breaking parameters are constrained at the 
gauge coupling unification point Mx- Instead, electroweak baryogenesis is not viable 
in a minimal SUGRA scenario of the MSSM. It requires an unconventionally large 
hierarchy between the left-handed and right-handed top squarks [75] , which is difficult 
to obtain within the framework of minimal SUGRA. In addition, the CP-odd soft 
phases required for successful electroweak baryogenesis face severe constraints from 
the absence of observable 2-loop contributions to the electron and neutron electric 
dipole moments [76]. 

(b) As has been discussed in Section 6, if one assumes that the neutrino- Yukawa cou- 
plings h\j have a certain hierarchical structure controlled by the approximate break- 
ing of global flavour symmetries, the model can have further testable implications for 
low-energy observables of lepton flavour and/or number violation, e.g. decay, 
p — > 67, p — > eee, — > e conversion in nuclei etc. In addition, electroweak-scale 
heavy Majorana neutrinos may be copiously produced at high-energy colliders, such 
as the LHC, the ILC and e~7 colliders, whose decays give rise to distinctive signa- 
tures of lepton-number violation which are usually manifested by like-sign dileptons 
accompanied by hadron jets. 

(c) Since successful inflation requires small couplings, i.e. k,, X, p < 10^^, the inflaton 
fleld decouples effectively from the low-energy spectrum and the Higgs-sector of the 
model becomes identical to the one of the MSSM. In spite of the aforementioned 
decoupling of the inflaton, however, the Fp-term hybrid model could still point to- 
wards particular benchmark scenarios of the MSSM. For example, ii X ^ k, the 
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FD-term hybrid model may explain the origin of a possible large value of the fx- 
parameter. Specifically, if A = 2k,, = —as = 2A'/susy! one gets from (2.4) the 
hierarchy /x ^ 4Msusy, where Msusy is a common soft SUSY-breaking scale of all 
scalar fermion fields in the model. If one additionally requires At = Ai, = 2Msusy; 
the low-energy fimit of the Fo-term hybrid model becomes identical to the so-called 
CPX benchmark scenario [77] describing maximal CP violation in the MSSM Higgs 
sector at low and moderate values of tan (3. In the CPX scenario, the lightest neutral 
Higgs boson weighing less than 60 GeV might have escaped detection at LEP. There 
have been several strategies to unravel the existence of such a fight CP- violating Higgs 
boson [78]. 

(d) The possible CDM scenario with the right-handed sncutrinos as LSPs requires large 
A and p couplings that could make Higgs bosons decay invisibly, e.g. H — > Ni^sp ^lsp- 
Also, right-handed sneutrinos could be present in the cascade decays of the heavier 
supersymmetric particles. The collider phenomenology of such a CDM scenario lies 
beyond the scope of the present article. 

The Fo-term hybrid model studied in this paper should be regarded as a first attempt 
towards the formulation of a minimal Particle-Physics and Cosmology Standard Model, 
which does not involve excessive fine-tuning in the fundamental parameters of the theory. 
As we outlined above, it might be possible to test the validity of our model by a number 
of laboratory experiments and further substantiate it by future astronomical observations. 
The Fo-term hybrid model is not plagued with the usual gauge-hierarchy problem of non- 
supersymmetric theories and can, in principle, be embedded within an E(6) or E(7) GUT, 
within the framework of SUGRA where SUSY is softly broken at the TeV scale. In the 
same vein, we note that a possible natural solution to the famous cosmological constant 
problem will shed valuable light on the model-building aspects of cosmologically viable 
models. It will also open up new avenues in quantitatively addressing the major energy- 
density component of today's Universe, the so-called Dark Energy. We hope that all these 
insights, along with new observational and experimental data, will help us to improve 
further our present bottom-up approach to formulating a more complete minimal model of 
particle physics and cosmology. 
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A Mechanisms of Explicit £)-Parity Breaking 



Here we will present mechanisms for explicitly breaking D-parity within the SUGRA frame- 
work, pointing out their possible implications for the decay rates of the (^-sector particles. 
We separately discuss the breaking of D-parity for an Abehan U(l)x and an non-Abelian 
SU(2)x waterfall-gauge sector. 

A.l Z)-Parity Breaking in the U(l)x Waterfall- Gauge Sector 

As already discussed in Section 4.1, the simplest way of breaking D-parity is to add a 
subdominant bare FI D-term Cpi to the Lagrangian [cf. (4.5)]. As was shown in [13], 
however, even if such a term were absent from the tree-level Lagrangian, it could still 
be generated by quantum-mechanical effects in an effective manner, after integrating out 
Planck-scale degrees of freedom. It should be stressed here that the radiative generation of 
an effective FI D-term occurs only after the SSB of the U(l)x gauge symmetry. 

To elucidate this point, let us consider a simple extension of the Fo-term hybrid 
model, which includes a pair of superfields Xi^2 of opposite U(l)x charge, i.e. QiX,) — 
—Q{Xi) = Q{Xi) = —Q{X2) ~ 1- In this case, the superpotential l^iw pertinent to the 
inflaton- waterfall sector may be extended as follows: 

= kS(x,X, - M') + im^X% + 6^4^ + ■■■ 
V / 2mpi 2mpi 

where the dots stand for subleading terms that multiply the leading operators by extra 
powers of [XiX,)^ /m^, with n > 1. These subleading operators arc irrelevant for our 
discussion here and can be ignored, within a perturbative framework of SUGRA. The 
leading operator form of the superpotential (A.l) may be reinforced by the R symmetry: 

S ^ e"^S, ii,2 ^ e^"/'Xi,2, {L, Q) ^ e'" (L, Q) , (A.2) 
with W e^°'W. As before, all remaining fields do not transform under the R symmetry.^ 
We will now show that a FI Z^-term, —^grupiD, will be generated if the ultrahcavy 
Planck-scale superfields Xi^2 are integrated out. As a starting point, we consider the U(l)jss: 
D-term Lagrangian 

Cn ^ ^D^ + - X*X2 - XX^ + Xx,) ■ (A.3) 

^Observe that the i?-syinmetry (A.2) allows for the subleading operator n' S{XiX2Y /rn^y This super- 
potential term can break the U(l)x gauge symmetry along the inflationary trajectory, thereby inflating 
away unwanted topological defects [25]. This scenario is known as shifted hybrid inflation. 
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Figure 11: Radiative generation of an effective FI D-term, — | mpj 

In order to integrate out the fields Xi^2, we need their mass spectrum in the post-infiationary 
era, where {Xi^2) = M and {Xi^2) = (5') = 0. In the weak basis X± = [Xi ± X2), this 
is approximately given by the Lagrangian 



->c^is ^ (^+, x*_) 



2mpy 




(A.4) 



The approximate mass eigenstates derived from (A.4) are 

X+ = X+ + seX., X_ = X_ - S0X+, (A.5) 

where sq (^1 — ^[)M'^/ (4^m|i) is a mixing angle which is typically much smaller than 1. 
In terms of the mass-eigenstates X±, the part of the Lagrangian (A. 3) hnear in the D-terms 
associated with the Planck-scale degrees of freedom reads: 

4" = -|/^(x;x_ + Tx, 



x\x_ + X*_X+ + 2s0 [x\x+ - x\xS) + 0(5^)] . (A.6) 



A FI £)-tadpole can only be generated from terms linear in se in the Lagrangian (A.6). 
This result should be expected on symmetry grounds, since terms linear in sg explicitly 
break the D-symmetry. The D-tadpole m|j, calculated from the one-loop graph of Fig, 11, 
is found to be 



Typically, one gets mpi/M < IQ-^, for M = GeV and ^1, < 10"^ 

For later use, it is interesting to outline a short-cut derivation of the result (A. 7), 
using the original weak basis of the fields, i.e. Xi 2 and Xi^2- We notice that, after the SSB 
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Figure 12: Diagram pertinent to a short-cut derivation of the effective FI D-term. 
of U(l)x, the F-terms of Xi^2 give rise to the D-odd operator, 

in the scalar potential of the extended Fo-term hybrid model. This operator induces, via 
the diagram shown in Fig. 12, an effective FI D-tadpole. Since the scalar fields Xi^2 are 
degenerate in mass to leading order, i.e. M-^^^ k, ^mpi, the graph in Fig. 12 is easily 
calculated using standard field-theoretic methods. It is logarithmically divergent, and in 
an effective cut-off theory it is given by (A. 7). We will use this short-cut approach below to 
calculate effective D-tadpoles in more elaborate extensions of the inflation-waterfall sector. 

The size of the FI D-term may be further suppressed, if the Planck-mass chiral 
superfields Xi^2 possess higher U(l)x charges. In general, one may assume that the U(l)x 
charges of Xi^2 are: Q{X2) = —Q{Xi) ~ n, where n > 1. In this case, the leading operator 
form of the inflaton- waterfall superpotential reads: 

Employing the short-cut method outlined above, it is straightforward to compute the loop- 
induced D-term, which is given by 

To obtain a small ratio mpi/M ~ 10^® with ~ 1, one would simply need n = 5, 6. 

Finally, we should remark that the loop-induced D-term does not lead to spontaneous 
breakdown of global super symmetry. 



A. 2 Z)-Parity Breaking in the SU(2)x Waterfall- Gauge Sector 

Here we outline two possible mechanisms for explicitly breaking the D-parities, Di and D2 
defined in (4.20) and (4.21), which govern the minimal inflaton-waterfall sector based on 
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an SU(2)x gauge group. 

The first mechanism utihzes a non-minimal Kahler waterfall-gauge sector, where the 
two -D-parities are broken by non-renormalizable operators. To be specific, a minimal 
Di,2-parity violating Kahler potential of the waterfall-gauge sector may be cast into the 
form: 



{Xj e-^s^^X*f k[ (Xle^s^^ir^Xs) (Xje^^^^^Xi) + H.c. 



+ «2 7: 9 + 



2 2 



4 {Xje^aVx j^^ ) (X^^e-^g^^ X*) + H.c. 
2m2 



^pi 

where the ellipses denote possible higher-order non-renormalizable operators. The couplings 
Ki^2 a-re real, whereas k'i 2 can in general be complex. Moreover, the difference = Ki — K2 
signifies I?i-parity violation, whilst parameter of D2-parity violation. 

Hence, non-zero values of the parameters k_ and k/_ will give rise to Di- and D2-parity 
violation in the waterfall-gauge Kahler potential i^^wF- Notice that, as far as Di-parity 
violation is concerned, the present mechanism applies to the Abelian case as well. 

There are several sources of D-parity violation contained in -ft'wF- More explicitly, 
D-paritj violation will first originate from the terms linear in D"", where are the aux- 
iliary SU(2)x components of the gauge-vector superfield Vx- In fact, these are the lowest 
dimensional Di,2-odd operators that emerge after the SSB of the SU(2)x and are given by 
the effective Lagrangian 

^D«-tad = 9_M^ ( ReK'_D^ - W_L>2 + k_ D') . (A.12) 
^ 2 mf,i V / ' ' 

These effective FI D-terms can be included in the Lagrangian by adding the constants 
2 ("^Fi)^ to the on-shell constrained terms, according to the scheme: — > D"-|-|(mpj)^, 
where 



[m 



Reft;'_, (mliY = 5- Im«;'_ , (mLf = H- . (A.13) 

mpj mpi 



One may obtain a fair estimate of the g^-sector particle decay rates, using the formula (4.18) 
and identifying mpi with mpj. In this way, we obtain 

r[-R_, -/+, ^ [kI, Re'{K'_), lm\K'_)] — . (A.14) 

IZOTT mpj 

In addition to the effective £)-tadpoles, higher-dimensional operators will also break the 
D-parities and so render the g^-sector particles unstable. For example, after expanding 
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the supcrficlds Xi^2 about their VEVs in the Kahler potential (A. 11), we find the non- 
renormalizable D-parity violating interactions described by the Lagrangian 



C. = 



non-ren — „ 9 



^ 2 + -7=-^ + H.c.)|a/X+|^ (A.15) 

^?7T.pi 4v2mpj V / 

With the aid of (A.15), an order of magnitude estimate of the g'-sector particle decay rates 
gives: Tg ~ |k'_P) g^M^ /m%i. These are of comparable order to the ones obtained 
earlier in (A. 14). For a typical inflationary scale, M — 10^^ GeV (with ~ 1), we find the 
decay width Tg ~ (k^, |«;'_|^) 10^ GeV. The latter should be compared with the bounds: 
3.8 X 10~^^ GeV < Fg < 4.3 GeV, corresponding to an upper and lower hmit on the second 
reheat temperature Tg of cosmological interest: 0.3 TeV ^ Tg ^ 10^ GeV. Consequently, 
values ranging from 10~^ to 10~^ for the couplings k_ and/or k'_ are required for suc- 
cessful coupled reheating. The lower end values of order 10~^ may possibly be seen as a 
strong tuning of the parameters. One way to explain the smallness of these parameters 
is to contemplate Kahler manifolds that break the D-parities by even higher-order non- 
renormalizable operators, e.g. of order ~ 1/mpj. In this case, the couplings k_ and k'_ will 
be multiplied by extra factors of M'^/rrip^ ~ 10^"^, so these couplings can have values of 
order 1 and still predict a second reheat temperature Tg ~ 0.3 TeV. 

We now describe a second mechanism of D-parity violation which might be useful 
to obtain small D-parity violating interactions. Let us therefore assume that the Kahler 
potential respects the D-parities. In this case, we may invoke a mechanism very analogous 
to the Abelian case. We extend the field content of the infiaton-waterfall sector by adding a 
pair of Planck-mass chiral superfields Xi and X2, which belong to the anti-fundamental and 
fundamental representations of SU(2)x, respectively. As in the U(l)x case, the superheavy 
superfields Xi 2 are charged under the continuous i?-symmetry given in (A. 2). With this 
restriction, the leading operator form of the infiaton-waterfall superpotential is given by 



iw 



kS^X.^X^ - M') + im^Sl% + e, 



mpi 



^ „ iX^VX2) jX.VX^) ^ ^ iX,VX2) (X2^^2) 
+ U2 + (,1 

mpi mpi 
mpi 

where the dots stand for additional operators that turn out to be irrelevant for the genera- 
tion of effective D"-tadpoles, and especially for those related to the D^- and D^-terms. The 
presence of these operators is only important to lift an accidental global U(l)j^ symmetry 
that governs this restricted part of the superpotential Wrw under consideration. Here, 
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T-1 



T-2 



✓ "X'^v ^'X'"^. ^'X^- 

/ \ / \ / \ 

-^1(2)^ f^l(2) ^1(2)^1 *^1(2) ^1(2)^ 1^^1(2) 



Figure 13: Diagrams responsible for the generation of effective D^'^'"^ -tadpoles for the SU(2)x 
case, in the single insertion approximation of the D-odd operators The subscripts in 

parentheses label the {Di,D2) parities of the respective operator. 



all non-renormalizable couplings ^1,2 and Ci,2 can in general be complex. Extending the 
notion of Di^2 parities to the Planck-mass superfields Xi^2, we observe that the operators 
related to the couplings k, ^ and 6*1,2 are even under Di and D2, whereas those related to 
the couplings Ci,2 are D2-odd. Moreover, the superpotential operators proportional to the 
couplings C+(-) = Ci + (~) C2 are Di-even (Di-odd). 

To calculate the effective D^'^'^-tadpoles after the SSB of the SU(2)x gauge group, 
we use the short-cut approach described above in Section A.l. Thus, the F-terms of Xi^2 
give rise to the following D-odd contributions to the scalar potential: 



0*2 C- 



2mpi 



-T 



(1-2)] 
- (1 - 2) 



+ H.c. 



(A.17) 



-772 



2 mil 



(xl(Xn) {x^^iT'{X2)) + il^2)_ 



2mpj 



x\ir'{X;)) fxf (Xi) ) + (1 ^ 2) + H.c. , (A.18) 



M2 



PI 



- a 



2) 



(xUT^x;))[{xl)zr^x,) + (x[{xi)) [{X^)X,) 

(A.19) 



where the subscripts in parentheses indicate the {Di,D2) parities of the above operators. 
Note that possible D-odd operators proportional to ^61^2 and have not been displayed, 
since they do not contribute to the generation of effective D"-tadpoles. To be specific, 
the effective D-tadpoles are induced radiatively via the graphs shown in Fig. 13, once the 
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operators JT^'^'-^ are individually contracted with the D-term operator T>^'^ related to the 
Xi^2 fields: 




These loop-induced effective FI D-terms can be included in the effective Lagrangian by 
shifting the on-shell constrained terms by constants, according to the above described 



where 9± — 9i:k 62. It can be estimated from (A.21) that for values ^± , C± ~ 10~^, one 
gets rripi'^/M < 10~^, leading to a low second reheat temperature Tg, below 1 TeV. In this 
context, one should notice that the size of the effective D-tadpoles is very sensitive to the 
cut-off scale, which we have chosen here to be the reduced Planck mass mpi. For instance, 
if a cut-off larger by one order of magnitude were adopted, then values of order 10~^ for the 
non-renormalizable superpotential couplings would be sufficient to generate the effective 
D"-tadpoles at the required size. 

The violation of D-parities will also affect the particle spectrum of the SU(2)x 
inflaton- waterfall sector. This will depend on the particular choice of the non-renormalizable 
part of the superpotential and Kahler potential. Our intention is not to pursue this issue 
any further here, by putting forward a specific non-minimal SUGRA scenario. Instead, our 
goal has been to explicitly demonstrate the existence of at least two mechanisms, which 
utilize the non-renormalizable part of the Kahler potential or superpotential to break the 
D-parities at the required order of magnitude. 



scheme: D"" —>■ -\- |(mpj)^. In this scheme, the mass parameters (rripj)^ are found to be 




(A.21) 
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